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Special Issue on High Polymers 


The following six papers were presented at the meeting of the Division of High Polymer Physics of the American 
Physical Society at Columbia University, January 30-February 1, 1947. 


Measurements of the Refractive Index of Films* 


FRED W. BILLMEYER, JR. 


E. I. du Pont de Nemours and Company, Plastics Department, Arlington, New Jersey 


The refractive index of a film of transparent or translucent material may be determined by 
immersing the film in a suitable liquid mixture and observing the intensity of light reflected 
from the surface of the specimen. The composition of the immersion liquid is varied, and a 
record is made of its refractive index (measured on a refractometer) and of the reflected light 
intensity. The refractive index of the film is the same as that of the immersion liquid giving the 
lowest reflected intensity. This has been confirmed by experiments with a glass plate, whose 
refractive index was known. According to Fresnel's laws of reflection, the experimental data 
should fit a parabola with its minimum at the refractive index of the film. On fitting the data 
by the method of least squares, the minimum is obtained with a standard deviation of the 
order of 0.002 refractive index units. The method has been applied to several plastic materials. 


INTRODUCTION 


METHOD has been developed for meas- 
uring the refractive index of a solid material. 
Originally used with thin films of translucent 
plastics, it may be applied as well to thicker 
specimens, of more highly absorbing or colored 


sample is the same as that of the immersion liquid 

which gives the lowest intensity of reflected light. 
THEORY 

When a material of refractive index mo, im- 

mersed in a liquid of refractive index n not too far 


material. Mentioned by R. W. Wood! and used to 
some extent in the investigation of anomalous 
dispersion phenomena, the procedure consists in 
immersing the specimen in a suitable mixture of 
liquids and observing the intensity of light re- 
flected from its surface. The composition of the 
immersion liquid is varied, and a record is made 
of its refractive indices (measured by a re- 
fractometer) and of the corresponding intensities 
of the reflected light. The refractive index of the 


* Presented at the meeting of the Division of High 
Polymer Physics, American Physical Society, January 30 
February 1, 1947. 

'R. W. Wood, Physical Optics (The Macmillan Com- 
pany, New York, 1934), third edition, pp. 411-12, 509. 


different from mo, is illuminated with unpolarized 
light striking its surface at an angle 7 (see Fig. 1), 
the fraction of the light reflected is given by 
Fresnel’s Law as? 


Fic. 1 








2 The coefficient } in these equations becomes nearly 1 
when reflections from the back face of a transparent sample 
are included. For these experiments its actual value is 
immaterial. 
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I ifsin?(a—r) tan*(¢—r) 
| L 


Iy 2ksin*(i+r) tan*(i+r) 
where /,) and 7 are the initial and reflected light 
intensities, respectively, and «@ and r are the 
angles of incidence and refraction. The angle r is 
given by Snell's Law: 


n sint = No sinr. (2) 


If 7 is chosen as 45°, a convenient angle experi- 
mentally, and m is nearly equal to mo, r becomes 
nearly 45°, and Eq. (1) reduces to 


I/Ig=}(tan*?(t—r) —tan*(t—r)). (3) 


Use of Eq. (2) then gives, for i—r small, 


| n—Ny\* 
I n=1(—) +--- (4) 
No 


which is a parabola with its minimum at n=». 
EXPERIMENTAL 


Apparatus 


A diagram of the essentials of the apparatus is 
shown in Fig. 2. 
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In the actual experiments two sets of apparatus 
were used. One consisted of a Bausch and Lomb 
projection microscope illuminator as light source 
and a Photovolt photometer, model 500, as 
detector. The tungsten lamp was run from a 
voltage-stabilizer, and a Wratten 62 filter was 
used to approximate mercury green light. 

Later experiments were made with the light- 
scattering apparatus built by the author (to be 
described elsewhere). A stabilized incandescent 
lamp was used in conjunction with Wratten 
filters 11 and 22 to approximate sodium yellow 
light. Occasionally the filters were omitted. 

In each case the sample was a thin film held 
between two pieces of polythene approximately 
0.080 inch thick, arranged so as to fit tightly 
along the diagonal of a 4-centimeter square glass 
cell. The polythene was cut away so as to allow 
the light to hit the film and to allow free circula- 
tion of the immersion liquid. Polythene is 
especially suitable because of its general in- 
solubility. 

The immersion liquids must be miscible with 
one another and must not dissolve, swell, or 
otherwise affect the sample. Some liquid pairs 
which have been used are benzyl alcohol and 
carbon tetrachloride for the range in refractive 
index around 1.5, and carbon tetrachloride and 
benzene, or carbon tetrachloride arid ethanol, for 
lower indices. 

Measurements were carried out by starting 
with a liquid mixture about 0.01 refractive index 
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unit away from the index of the film (chosen by 
experience or in a preliminary experiment) and 
adding small increments of a liquid with index 
beyond that of the sample. For each liquid mix- 
ture used, the refractive index and the intensity 
of reflected light were measured. 

The measurements of refractive index were 
carried out either on a Bausch and Lomb ‘“‘Pre- 
cision” refractometer using mercury green light, 
or on a Bausch and Lomb Abbé refractometer, 
which gives a reading for sodium light. All ex- 
periments were carried out at room temperature 
(approximately 25°C) and the refractometers 
were checked against standard glass plates 
occasionally. 

Glass Plate Experiment 

To test the accuracy of the procedure, the 
refractive index of a microscope slide was meas- 
ured both by the method described here and 


directly by the Abbé refractometer, with the 
following result: 


by refractometer, 
n=1.5171 
by light reflection, 


n= 1.5156+0.0036 (95 percent reliability). 


The two values agree within one standard 
deviation. The data are presented in Fig. 3, with 
the refractometer value indicated by an arrow. 
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TABLE I. Refractive indices of films. 








Material Index Color of light Notes 
polythene 1.515 mercury green PM-1, lot 720 
1.519 mercury green PM-1, lot 859 





tetrafluoroethylene 1.378 white ‘*Tefion,"’ density 2.18 
1.376 sodium yellow “Teflon,’’ density 2.12 


hexamethylene 1.532 sodium yellow Type FM-3 Nylon 
sebacamide 


COMPUTATIONS 


The refractive index of the test specimen may 
be obtained from the experimental data either by 
inspection or from the equation of the parabola 
fitted to the observed points by the method of 
least squares. 

In the latter case, the reflected intensity is 
represented by the equation 


IT=a+bn+cn’? (5) 


and the position of the minimum is 


no=-— (6) 
2c 


In the experiments carried out on films of 
plastics, the standard deviation in the value of mo 
was computed directly from the experimental 
data for the case of the glass plate. For each of 
the other cases the standard deviation was esti- 
mated from this value and the relative magni- 
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tudes of the “‘chi-square’’ index of precision® for 
glass and for the case in question. 

Whenever the refractive index as determined 
by a visual inspection of the data was compared 
with that computed by least squares, the differ- 
ence between the two was less than twice the 
standard deviation of the computed value. 

Some experimental results are presented in the 
following table, and a few of the computed curves 
are shown in Figs. 4, 5, and 6. Included on the 
graphs are the estimated and computed indices of 
refraction of the respective films, and a scale of 
standard deviation, ¢. The intensity is represented 
in arbitrary units, with the zero of intensity con- 
siderably below the axis. 

The intensity at the minimum will depend 
upon the transparency of the film. For a very 


3 See, for example, R. A. Fisher, Statistical Methods for 
Research Workers (Oliver and Boyd, 1941), eighth edition. 


transparent material such as glass, it will be 
limited chiefly by the scattered and stray light in 
the apparatus. The value of mo is not influenced 
by the magnitude of this residual intensity. 


CONCLUSION 


A method similar to the mineralogical im- 
mersion technique for measuring the refractive 
index of a solid has been presented. The test 
specimen may be a thick or thin piece with a 
reasonably plane face, or a film. The material 
may be clear or translucent, and may also be 
colored. The optical quality of its surface need 
only be fair. For these reasons this technique has 
proved successful in’ several cases where other 
methods, such as the determination of Brewster’s 
angle, have failed. 

Some experimental results are given in the 
accompanying table and figures. 





A Sound Velocity Method for Determination of Molecular Weight of Liquid Polymers{ 


ALFRED WEISSLER, JAMES W. FITZGERALD,* AND IRVING RESNICK** 
Naval Research Laboratory, Washington, D. C. 


A new method is presented for determining number 
average molecular weights of liquid polymers, by means of 
easily performed measurements. The molecular weight is 
an explicit function 


HE molecular weight of a polymer is so 

fundamental a property, and the methods 
available for its determination are so few in 
number, that any new technique warrants a 
certain interest. This paper presents such a new 
method, suitable for liquids, which is based on 
the relatively simple measurements of sound 
velocity, density, and refractive index. 


I. BASIS OF THE METHOD 


Sound velocity in liquids became useful for 
the investigation of molecular structure after 


+ Presented at the meeting of the Division of High 
Polymer Physics, American Physical Society, January 30- 
February 1, 1947. 

* Now at Engineering Research Associates, Washington, 


D. ¢€ 
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of the sound velocity v, density d, refractive index n, and 
the two empirical constants A and B. Accuracy of about 2 
parts per hundred was attained for the lower polyethylene 
glycols. The method seems suitable for molecular weights 
up to several thousand. Sound velocities were measured 
by an acoustic interferometer, at one megacycle frequency. 


the discovery by M. R. Rao! of the empirical 
constant 


R= Mov!/d (1) 


where M is the molecular weight, v and d are the 
sound velocity and density measured at the same 
temperature, and R may be called the molar 
sound velocity. For each pure liquid (except 
water) the molar sound velocity is a constant 
independent of temperature; therefore, density 
and the cube root of velocity vary in such a way 
that their ratio is fixed. 





** Now at Johns Hopkins University, Baltimore, 
Maryland. 

1M. R. Rao, Ind. J. Phys. 14, 109 (1940); J. Chem. Phys. 
9, 682 (1941). 
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TABLE I. Physical properties of polyethylene glycols at 30°C. 























Sound Boiling 

Glycol Mol. wt. Density velocity nD Viscosity point 
Ethylene 62.1 1.1068 g/cc 1642.9 m/s 1.4295 13.95 c.p 69°/1.2 mm 
Diethylene 106.1 1.1094 1567.7 1.4438 22.55 105°/2.0 
Triethylene 150.1 1.1158 1593.3 1.4531 29.96 123°/1.2 
Tetraethylene 194.2 1.1159 1580.1 1.4567 35.16 138°/0.5 
Pentaethylene 238.2 1.1165 1580.0 1.4593 42.37 176°/0.5 
Poly 200 200+ 10 1.1168 1592.2 1.4572 39.71 
Poly 300 300+15 1.1176 1578.0 1.4618 59.02 
Poly 400 400+20 1.1182 1576.0 1.4636 74.70 
Poly 600 600 +30 1.1183 1570.2 1.4653 104.62 

Although the velocity of sound varies irregu- molecular weight 

larly within a given homologous series, such as vv Bd 7 
the alcohols or ketones, the molar sound velocity aalie (5) 


was found to be an additive property of the atoms 
(or bonds) in the molecule. From the published 
tables of R increments for atoms of carbon, 
hydrogen, oxygen, etc., the molar sound velocity 
may be calculated with good accuracy, just as in 
the familiar case of molar refraction. 

Further, Lagemann and Dunbar? showed that, 





n*—1 
vi-A 
Ga 
in terms of the two empirical constants and the 
observed sound velocity, density, and refractive 
index. Similar expressions may be obtained from 
other pairs of molar constants; for example: 








within a homologous series, linear relationships B'd 
exist between any two of the following molar M= nem (6) 
: ‘ v!—A’(2.9+log logn) 
constants: (a) molar sound velocity, (b) molar _ 
refraction . B’d 
Myn?-1 hata (7) 
Pe am : pa, 
d \n?+2 re —A’'(2.9+log logn) 
n*+2 


where 1 is index of refraction, (c) molar viscosity 
as defined by Souders?® 


where A’, A’, B’, and B” are the analogous 
slopes and intercepts of the other linear relations. 


Vv These results should’ be applicable not only to 
I= (2.9+log logn)— (3) & homologous series, in which the difference be- 
t tween two successive members is —CHe—, but 


where 7 is viscosity in millipoises, (d) molar 
rotation, (e) parachor, (f) 
volume, and (g) van der Waals’ b. 

An example of such a relationship is the one 
involving refraction and sound velocity 


magnetic critical 


R=AN+B (4) 


where A and B are, respectively, the slope and 
intercept. Upon substitution of Eqs. (1) and (2), 
this immediately yields an expression for the 


2 R. T. Lagemann and W. S. Dunbar, J. Phys. Chem, 49, 
428 (1945). 
3M. Souders, J. Am. Chem. Soc. 60, 154 (1938). 
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also to a polymer homologous series, where the 
difference is any monomer structural unit. 


II. EXPERIMENTAL 


It was therefore decided to test the validity of 
Eq. (5), and incidentally of (6) and (7), by using 
it to determine the molecular weights of several 
pure members and commercial mixtures of the 
polyethylene glycol series. These polymers, of 
type formula HO—(—CH,.—CH.—O-),—H, 
have been studied intensively in the past.*® The 

4H. Staudinger, Der Aufbau der hochmolekiilaren organ- 
ischen Verbindungen, (Julius Springer, Berlin, 1932). 


5 R. Fordyce and H. Hibbert, J. Am. Chem. Soc. 61, 
1910, 1912 (1939). 
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Fic. 1, Circuit design for interferometer. 


chemical individuals (the monomer, dimer, 
trimer, tetramer, and pentamer) were purified by 
fractional distillation through a 40 cm Widmer 
column, at about 1 mm pressure. The mixtures 
were those of average molecular weight 200, 300, 
400, and 600, as manufactured by Carbide and 
Carbon Chemicals Corporation. 

Several physical properties of these nine liquids 
were measured at 30°+0.1°: refractive index for 
sodium light by the Abbé refractometer, density 
by a 25 ml pycnometer, sound velocity by the 
acoustic interferometer described below, and 
viscosity by the Ostwald-Fenske viscometer. 
These properties are listed in Table I. 

The acoustic interferometer® was of con- 
ventional type,’~® energized by a one megacycle 
crystal-controlled oscillator. Figure 1 shows the 
circuit used. Figure 2 is a sectional view of the 
brass interferometer cell, containing the quartz 
crystal source and the movable reflecting plate, 
immersed in the thermoregulated oil bath. About 
300 ml of liquid are required for a measurement. 
As the reflector is moved vertically through the 
liquid by the micrometer head, cyclical variations 
in the standing wave pattern recur at distances 
of integral half-wave-lengths. These cause corre- 
sponding cycles in the current through the 
interferometer crystal, which are observed on a 
suitable microammeter. The micrometer-head 
travel required for an interval of two current 


® Designed by one of us (J. W. F.) and Mr. Burton G. 
Hurdle. 

7G. W. Pierce, Proc. Am. Acad. 60, 269 (1925). 

8 J. C. Hubbard and A. L. Loomis, Phil. Mag. 5, 1177 
(1928). 

®E. Klein and W. D. Hershberger, Phys. Rev. 37, 760 
(1931). 
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Fic. 2. Interferometer in thermoregulated oil bath. The 
inner tank shown has since been removed. 


maxima, then, is equal to the wave-length; and 
the velocity of sound is, of course, the product of 
this wave-length and the frequency. In practice, 
one measures the travel for twenty maxima, 
which makes possible a precision of a few 
hundredths of a pefcent. 


III. DISCUSSION 


From the physical properties of the nine liquids, 
there were computed the molar refractions, 
viscosities, and sound velocities. These observed 
molar constants agree reasonably well (Table I1) 
with the values calculated!’ from the molecular 
structures, using the published tables**'° of the 
bond increments for each type of constant. In the 
case of molar sound velocity, the agreement is 
better than 1 percent, which is less than the 
dispersion in values among the original com- 


TaBLeE II. Molar constants of polyethylene glycols. 


Molar sound velocity Molar retraction Molar viscosity 


Ob- Calcu- Ob- Calcu- Ob Calcu- 

Glycol served lated served lated served lated 
Ethylene 661.7 652.0 14.47 14.49 181.2 176.0 
Diethylene 1111.2 1106.1 25.40 25.39 312.9 316.9 
Triethylene 1571.6 1560.1 36.37 36.31 443.2 457.8 
Tetraethylene 2027.1 2014.2 47.36 47.28 575.3 598.7 
Pentaethylene 2485.1 2468.2 58.36 58.27 708.3 739.6 
Poly 200 2091 2075 48.79 48.87 593.6 617.5 
Poly 300 3125 3105 73.71 73.56 897.3 937.5 
Poly 400 4163 4135 98.63 98.48 1201.3 1257.5 
Poly 600 6236 6195 148.40 148.17 1813.4 1897.5 


10 Molar refraction was calculated by the method of K. 
Fajans, Physical Methods of Organic Chemistry, edited by A 
Weissberger (Interscience Publishers, New York, 1946), 
Vol. I, pp. 672, 679. 
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Fic. 3. Linear relation between molar sound velocity and 
molar refraction. 


pounds used for constructing the table of incre- 
ments. Compared to values previously reported" 
for 20°, the molar refractions show the customary 
increase of 0.01—0.02 percent per degree, at- 
tributed by Fajans!° to loosening of the electron 
shell. The agreement is poorest for the molar 
viscosities. 

Next, the existence of a linear relation between 
molar sound velocity and molar refraction was 
verified for the five pure compounds. The points 
lie on the straight line of Fig. 3, of which the 
slope A is 41.59 and the intercept B is 57.7. The 
slope may also be calculated from the tables of 
bond increments: for each —CH2:—CH:—O— 


TABLE III. Determination of molecular weight by 
sound velocity. 


v= t= ; 
v!— 41.59 
_ Ga) 
Mole. Mole. Mole. 
wt. by wt. wt. 
Theoretical sound number wt. 
Glycol Mole. wt. velocity Deviation avg avg. 
Ethylene 62.1 59.9 —3.5% 60.8 62 
Diethylene 106.1 111.5 5.1 101.6 109 
Triethylene 150.1 147.3 -—1.9 148.5 155 
Tetraethylene 194.2 195.7 —0.8 195.4 190 
Pentaethylene 238.2 237.7 0.2 234.5 238 
Poly 200 200+10 186.2 —2.0+ 183 219 
Poly 300 300+15 302.7 0.0+ 315 359 
Poly 400 400+20 381.8 0.0+ 384 480 
Poly 600 600+30 546.9 -—4.0+ 561 720 


11 A, F. Gallaugher and H. Hibbert, J. Am. Chem. Soc. 
58, 813 (1936). 


VOLUME 18, MAY, 1947 



































3000, 60 
E ? 
. z 
8 S 
w@ 200 S 
rm 409 
2 th 
rod c 
ed « 
& $ 
3 100 20$ 
4 
200 400 600 B00 «(OO 


MOLAR VISCOSITY 


Fic. 4. Verification of other linear relations. 


structural unit, the R increment is 454 and the N 
increment is 10.88. The quotient of these two, 
which should be the slope, is 41.73, so that the 
agreement is a few parts per thousand. 

As expected, linear relations were found (Fig. 4) 
for the two other pairs of molar constants, sound 
velocity and viscosity, and refraction and vis- 
cosity. For the former pair, slope A’ = 3.466 and 
intercept B’=31.9; for the latter, A’’ =0.08337 
and B’ = —0.640. The agreement between ob- 
served and calculated slopes is not as good in 
these cases. 

Next came the crux of the entire work. For 
each of the nine liquids, the molecular weight was 
calculated as indicated by Eq. (5), from the 
density, sound velocity, index of refraction, A, 
and B. The results show an average accuracy of 
about 2 percent (Table III) on the basis of the 
theoretical molecular weights, which compares 
well with that of other methods. 

For purposes of comparison, the number aver- 
age molecular weight was determined by freezing- 
point depression in water as solvent.'®? Meas- 
urements were made at concentrations of 
approximately 2 percent and 4 percent, and the 
values extrapolated to infinite dilution. The pre- 
cision is estimated at about 3 parts per hundred. 
It is apparent from Table III that, within ex- 
perimental error, the sound velocity molecular 
weight is the same as the number average 


12 Using water as solvent avoided the anomalous results 
encountered in ethylene dibromide solution by Gallaugher 
and Hibbert, reference 11, presumably because the glycols 
tend to be solvated, rather than associated, in aqueous 
solution. The apparent molecular weights obtained in the 
present work decreased slightly with increasing concentra- 
tion, which indicates that association of glycol molecules is 
negligible compared to solvation. 
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Fic. 5. Bulk viscosities of polyethylene glycols at 30°. 


molecular weight. In addition, a weight average 
molecular weight was determined from the log- 
log plot of bulk viscosity against molecular 
weight" (Fig. 5). As would be expected, for the 
heterogeneous mixtures the weight average 
results shown in Table III are considerably 
higher than the number average or sound 
velocity molecular weights. 

Attempts were also made to calculate the 
molecular weights from Eqs. (6) and (7), the 
analogous relations involving the other two pairs 
of properties. However, the errors are large, as 
demonstrated in Table IV. This may be due 
partly to the fact that the precision of measure- 
ment of viscosity is only about one-tenth as good 
as that of sound velocity or refractive index. 
Unfortunately, the errors increase for higher 
molecular weights because the calculations in- 
volve a decreasing difference between the two 
terms in the denominator of the expression for M. 


‘3 These lower polyethylene glycols appear not to conform 
to the relation logy = A+B M! suggested for bulk viscosities 
by P. J. Flory, J. Am. Chem. Soc. 62, 1057 (1940). 
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IV. SUMMARY 


A new method for determining molecular 
weights of polymeric liquids has been shown to 
yield results of about 2 percent accuracy in the 
case of the lower members of the polyethylene 
glycol series. This method, which gives the 
number average molecular weight, consists of 
calculating M from two empirical constants and 
the observed sound velocity, density, and re- 
fractive index. We intend to investigate extension 
of its applicability to solid polymers or solutions 


_of solids. . 
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Refraction Velocity Refraction 
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An Instrument for Measuring Particle Diameters and Constructing Histograms from 
Electron Micrographs* 


E. E. HANSON AND J. H. DANIEL 
Chemical and Physical Research Laboratories, The Firestone Tire and Rubber Company, Akron, Ohio 


An instrument for shortening the time and labor required in constructing particle size 
histograms from electron micrographs is described. Representative histograms obtained with 
the instrument for GR-S latices are shown. Linear, surface, and volume-equivalent diameters of 
non-spherical particles are defined, and variations in the instrument for plotting such diameters 


are proposed. 


INTRODUCTION 


N certain studies of the polymerization of 

synthetic rubbers, it is desirable to obtain 
particle-diameter distribution curves or histo- 
grams for the rubber latices. Such diameter 
histograms may be obtained from electron micro- 
scope photographs of the latex particles. 

The technique used in this laboratory for 
preparing latex specimens for the electron micro- 
scope! consists essentially of forming a film 
from a dilute water solution of polyviny! alcohol 
containing a small amount of the latex. 

The construction of particle diameter histo- 
grams is tedious and time consuming. In the 
usual procedure, the “longest”? and “‘shortest”’ 
diameters of the projected image (or of a photo- 
graphic enlargement) of a particle are measured 
with a rule, their arithmetic mean is computed, 
and finally the number of particles with mean 
diameters lying in chosen diameter intervals is 
plotted against the diameters. Multiplication by 
a scale factor, containing the reciprocal of the 
magnification, converts the diameter axis of the 


histogram to read directly in Angstrom units. 
To obtain a representative distribution, at least 
200 particles should be measured. 

The instrument to be described here was de- 
veloped to reduce the time and labor required in 
constructing such histograms. It will be called a 
“histogram computer.” 

In the development of the histogram com- 
puter, consideration was given to the question 
of whether the arithmetic or the geometric 
mean, or some other function of the “‘longest’”’ 
and ‘‘shortest’”’ diameters of a particle, should be 
plotted. With an unbroken film in the electron 
microscope, the latex particles usually appear 
approximately spherical in shape. However, when 
holes are present in the film, as frequently is the 
case, the particle contours are elliptical. The 
problem, then, is to compute the diameter of the 
equivalent spherical particle having either the 
same volume, or the same surface area, as the 
actual ellipsoidal particle. The two cases require, 
of course, calculation of different functions of 


the major and minor axes. In addition, the 


TABLE I. Comparison of arithmetic and geometric mean diameters with volume and surface equivalent diameters. 


Prolate Spheroid 
I 


Oblate Spheroid 


Arith. 
Mean % Error % Error % Error ™ Error 
a+b Geom. Differ- 

Ratio - ean ence Arith. Geom. Arith. Geom. Arith. Geom. Arith. Geom. 
a/b 2 (ab)4 % De Mean Mean Ds Mean Mean D Mean Mean Ds Mean Mean 
1.0 b b 0.0 b 0.0 0.0 b 0.0 0.0 b 0.0 0.0 b 0.0 0.0 
1.5 1.255 1.2256 1.0 1.1455 9.2 7.0 1.166 7.8 5.6 1.315 —4.6 —6.5 1.335 —6.0 —7.9 
2.0 1.50b 1.4146 2.9 1.2605 19 12.2 1.31b 14 7.9 1.5875 —5.5 —1i11 1.66b -—9.6 —15 
2.5 1.755 1.585 5.1 1.357b 29 16 1.44b 22 9.7 1.842) —12 —21 





a= Major Axis 
b=Minor Axis 


* Presented at the meeting 
30—February 1, 1947. 


of the Division of High Polymer Physics, American Physical Society, New 


—5.0 —14 


1.995 


D,» =Diameter of sphere of equal volume. 
D, =Diameter of sphere of equal surface. 


York, January 


1R. H. Kelsey and E. E. Hanson, J. App. Phys. 17, 675 (1946). 
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Fic. 1. The scanning mechanism. 


particles may be either prolate or oblate, or 
some shape in between. Without going through 
some time-consuming process, such as taking 
stereo-photographs, there is no way of deter- 
mining what volumetric shape a particle has. 

A comparison between the diameters of the 
volume-equivalent and the surface-equivalent 
spheres, and the arithmetic and geometric mean 
diameters for both prolate and oblate spheroids 
is given in Table I. In some cases the geometric 
mean is the closer approximation, and in others 
the arithmetic mean is the closer. It was there- 
fore decided to design the simplest of possible 
mechanisms, which gives the arithmetic mean. 
In general, the measurements should be restricted 
to particles having a major to minor axis ratio 
of less than two. 


APPARATUS 


The histogram computer consists essentially 
of a scanning mechanism for setting manually 
two pairs of mutually perpendicular lines on the 
outline of the particle, and a coupled recording 
mechanism for dropping a steel ball into a 
compartment corresponding to a range of diam- 
eters containing the arithmetic average of the 
two measured diameters of the particles. 

Figures 1 and 2 are photographs of the scanning 
and the recording mechanisms, respectively. The 
scanning mechanism consists of a steel base- 
plate with a right angle section cut out for view- 
ing the electron micrograph field, two mutually 
perpendicular screws (one labelled E) for trans- 
lating non-rotating nuts to which clear Lucite 
blocks are attached, two equal gear trains 
coupling the rotations of the screws to their 
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Fic. 2. The recording mechanism. 


respective Selsyn motors, and a sub-base plate 
of clear Lucite on the under side of which are 
two scribed, inked-in lines drawn parallel with, 
and perpendicular to the two screws. Each of 
the movable Lucite blocks has a scribed, inked-in 
line drawn on its top and bottom surface, in a 
vertical plane parallel to a line in the sub-base 
plate. This enables the operator to set the lines 
without parallax on the edge of the particle. 
Each of the two Selsyn motors is coupled to a 
corresponding motor on the recorder so as to 
transmit its rotation. 

In the recorder unit, the Selsyn motors, A 
and B, transmit their rotations additively through 
a differential gear, C, to a screw 10 inches long. 
The screw serves to position the ball-dispensing 
mechanism £. Directly below the ball-delivery 
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Fic. 3. The ball release mechanism. 
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tube, and parallel to the screw, is a thin rec- 
tangular Lucite box, F, divided into 80 com- 
partments for receiving the steel balls as they 
are dropped. 

The ball-dispensing unit consists of a }’’ 1.D. 
copper-tubing magazine, D, for holding a supply 
of 34"’ diameter steel balls, and an electromag- 
netically controlled ‘“‘valve’’ by which the oper- 
ator can cause the balls to be dropped one at a 
time. A cross-sectional diagram of the valve is 
shown in Fig. 3. The axis of the delivery-tube is 
sufficiently offset from the magazine, D, to 
prevent the balls from falling through. The flat 
plunger, P, normally held by the spring, S, in 
the position shown, is pierced by a hole slightly 
larger in diameter than the steel balls. When the 
circuit to the magnet, M, is closed, the plunger 
carries one ball over and discharges it into the 
delivery tube. The delivery tube projects be- 
tween two parallel guides to a position 75” 
above the opening in the recorder compartments. 
The parallel guides also serve the purpose of 
preventing any rotation of the dispensing unit 
nut. 

The recorder compartment box, F, was made 
by bolting together the two ends of two pieces 
of 6”x11" xi” Lucite blocks against .096” 
thick spacers. Strips of .010” thick shim stock 
were inserted into grooves which had _ been 
accurately milled into the inner surfaces of the 
plates on }”’ centers. The bottom of the box is a 
metal blade which can be removed at the end of 
a run to release the accumulated balls into a flat 
box below. 

End-play was taken out on the three screws 
by special thrust bearings. Back-lash was elimi- 
nated through the use of split-nuts engaging 
the screws. 

The Selsyn motors were designed to be oper- 
ated normally at 115 volts and 400 cycles. 
However, since no 400 cycle source was avail- 








Fic. 4. Appearance of the 
scanning lines set on a 
particle. 








VOLUME 18, MAY, 1947 





Fic. 5. Contact print of the recorder compartmented box 
showing distribution from which histogram of run 3, Fig. 6, 
was obtained. 


able, they were operated from the transformer, G, 
at 60 cycles and 24 volts. 

The gear train normally used gave an amplifi- 
cation factor of 18 from the scanning to the 
recorder units. Therefore, a total magnification 
of 35,300 diameters was required in the electron 
micrograph to give a value of 50 Angstroms to 
each }’’-long compartment interval. Other values 
could be given to the compartment interval by a 
suitable choice of gears and total magnification 
of the electron micrograph. 

In practice, a particle count can be made 
either from a photographic enlargement, or 
directly from an image of the electron micro- 
graph plate projected at the desired magnifica- 
tion on a ground-glass screen. The procedure is 
first to set the scanner unit on the zero position 
(scribed lines directly over each other), and the 
ball dispensing unit approximately at the zero 
position of the recorder box, and finally to shift 
the compartmented box slightly by hand until 
the balls fall equally often on both sides of the 
chosen zero wall. The scanning unit is then set 
successively on the particles to be measured, so 
that each particle is bracketed by the scribed 
lines as shown in Fig. 4. As each particle in turn 
is thus bracketed, the ball release switch is 
closed to record that particle. A counter, coupled 
to the switch, automatically totals the number of 
particles as the measurements proceed. 
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Fic. 6. Histograms of experimental latex. Each run made 
from same electron micrograph field. 


TYPICAL RESULTS 


Figure 5 is a direct contact print of the record 
given by the histogram computer of 335 particles 
of an experimental GR-S latex. The print was 
made by placing F-2 contact paper against the 
back of the recorder compartments, and illumi- 
nating the front of the box for 15 seconds with a 
60 watt Mazda bulb placed 7 feet away. In this 
case the electron microscope magnification was 
5600, and the photographic enlargement was 
6.3 X to give a total magnification of 35,300 x. 

In Fig. 6 are three histograms of an experi- 
mental GR-S latex. All were constructed from 
the same electron micrograph field, and are, 
therefore, indicative of the reproducibility of the 
instrument. In these measurements an interval 
of 100 Angstrom units was chosen. Since the 
“limit of resolution of the electron microscope was 
about 50 Angstroms for the specimens used, it 
was concluded that the sharp variations between 
50 Angstrom intervals, shown in Fig. 5, were not 
significant. 

Figure 7 is similarly a record of three successive 
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Fic. 7. Histograms of another experimental latex. Each run 
made from same electron micrograph field. 


histograms constructed from the same field of 
an electron micrograph of another GR-S type 
latex. 

In general, the reproducibility of the curves is 
considered satisfactory. 

It has been our experience that it is much less 
tedious to operate the computer than to measure 
particles with a rule, that the accuracy is in- 
creased, and that the time required to construct 
a histogram is reduced by about one-half. 


OTHER POSSIBLE MECHANISMS 


Construction of an electronic servo-system, to 
replace the Selsyns and differential in the histo- 
gram computer, was not justified in this labora- 
tory because of extra complications, but would 
give added features which might be desirable in 
some applications. Among such features may be 
listed the following: 


1. Distribution of the geometric as well as arithmetic 
means of major and minor axes. 

2. Distribution of surface area (approximate) or of 
volume (exact), assuming particle shape to be an ellipsoid 
of revolution. 

3. Flexibility of size of histogram increments. 


One such scheme is shown in Fig. 8. Circular 
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potentiometers, P; and P», have resistances linear 
with their respective angles of rotation, @; and 4.2, 
to the desired degree of accuracy. In use, 6; 
and 62 are set proportional to the major and 
minor axes, D; and Dz, respectively. Similar 
potentiometers, P; and P,, are mechanically 
inter-connected so that their angles of rotation 
remain always equal to the single angle @. This 
angle 6 is determined by the position of the servo- 
motor, which also runs the ball dispenser of the 
recorder (minus differential gear). In Fig. 8, the 
two voltages, V2 and V3, applied to the sum 
amplifier will always differ 180° in phase, so 
when their magnitudes are equal, the sum am- 
plifier output to the servo-motor will be zero. 
The direction of rotation of the servo-motor 
when V2 and V3 are not equal will depend on 
which voltage is the greater. The blocks labeled 
CF are cathode followers with amplification, 
Acr,.extremely constant and slightly less than 
unity. Setting Ri=R2=0 for the moment, the 
following relations will exist: 


Ve= V»Acr)6.= V AcrykeD, 
=Rk,D,C\VoAcrik2De, (1) 
V3= VsAcr.6= —Acr26*CqVo, (2) 


where 0,;=k,D;; 02=k2.D2; voltage across P; 
=(C, Vo; voltage across Py= —C,Vo. 

When the servo-motor has equalized V; and 
— V2, and the recorder ball dispenser comes to 
rest, the following relation will exist: 


A CF207C4 | 0 = k,D,C, ] "0A CFiR2Ds», 
Ci AcF, i (3) 
d= (ut -—) (D,D,)}. 


C4 A CF2 





Hence the position of the recorder ball dispenser 
is proportional to the geometric mean of the 
particle major and minor axes. 

Since Vo does not appear in Eq. (3), line volt- 
age regulation is unnecessary. Rheostat R; allows 
the increment of diameter for the histogram to 
be varied. The combination of R; and R»2 allows 
a shift of recorder zero for cases where no small 
particles are present. 

To obtain a histogram of approximate surface 
areas, a switching arrangement feeds V4 instead 
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Fic. 8. Schematic diagram of proposed servo-mechanism 
for histogram computer. 


of V3; to the sum amplifier, thus giving @ pro- 
portional to the product D,D». 

To obtain a histogram of the arithmetic mean 
of particle diameters, a switching arrangement 
removes CF, and CF; from the circuit, connects 
the ungrounded ends of P; and P2, and feeds 
V; and V2 separately, and V4 instead of V3, to 
the sum amplifier (for this arrangement k; must 
equal ke, as would normally be the case). 

To obtain a histogram of volumes, V4 instead 
of V3 is fed to the sum amplifier, P3 is connected 
mechanically and electrically to P2 in the same 
manner it was formerly connected to P,, and the 
series outputs of P:, Ps, and P3, are fed to the 
sum amplifier. In the case of a prolate spheroid, 
6; must be set proportional to the major axis; 
for an oblate spheroid, @; must be set propor- 
tional to the minor axis. Obviously, only a small 
range of particle sizes can be covered in volume 
histograms as compared to diameter histograms. 
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Torsion of a Rubber Cylinder* 


R. S. RIVLIN 
British Rubber Producers’ Research Association, Welwyn Garden City, Herts., England, and 
National Bureau of Standards, Washington, D. C., U.S. A. 


It has been predicted theoretically that, in general, a right-circular cylinder of incompressible, 
highly elastic material, which is isotropic in its undeformed state, cannot be held in a state of 
pure torsional deformation by means of a torsional couple alone. In addition, normal surface 
tractions must be exerted over the plane ends of the cylinder. These normal surface tractions 
depend on the amount of torsion and on position on the plane ends of the cylinder. Experi- 
ments are reported here in which this phenomenon is observed in a right-circular cylinder of 
pure gum compound. The dependence of the surface traction on amount of torsion and its 
distribution over the surface of the cylinder is studied by measuring the bulging of the rubber 
into small holes in a metal plate on one end of the cylinder. ' 


I. INTRODUCTION 


N a series of papers which is pending publi- 
cation elsewhere,' the mechanical properties 

of highly elastic materials, such as rubber, which 
can be assumed incompressible and isotropic in 
their undeformed state, are discussed in terms of 
a stored-energy function W. This stored-energy 
function is specified in terms of the invariants of 
strain. It is shown that, once W is specified in 
this manner, the forces, which must be applied 
in order to support any given state of deforma- 
tion, can be calculated. This calculation has been 
carried out for the pure torsional deformation of 
a right circular cylinder, for an assumed stored- 
energy function of simplest possible type and 
subsequently for any general form of stored- 
energy function. 

The theory predicts that in order to produce a 
pure torsion in the cylinder, a torsional couple 
must be exerted on the plane ends of the cylinder, 
as in the classical theory of elasticity for small 
deformations. In addition, normal tractions must 
be exerted over these surfaces. The distribution 
of these normal tractions over the ends of the 
cylinder and their dependence on the amount of 
torsion have been calculated for any specified 
form of the stored-energy function. 

In this paper experiments are reported in 
which this distribution and dependence on 
amount of torsion is measured for a cylinder of 
pure gum compound, on the assumption that the 

* Presented at the meeting of the Division of High 
Polymer Physics, American Physical Society, January 30- 
February 1, 1947. 


'R. S. Rivlin, Phil. Trans. Roy. Soc. A (pending publi- 
cation). 
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surface traction is measurable by the bulging of 
the rubber into small holes in a metal plate on 
one end of the cylinder. The results obtained are 
compared with those derived theoretically for 
two simple forms of the assumed stored-energy 
function. Agreement, within the accuracy of the 
experiments, is obtained if the stored-energy 
function has a form equivalent to that put 


forward for rubber by Mooney.? This form in- 
volves two physical constants for the material. 


II, SPECIFICATION OF THE STORED- 
ENERGY FUNCTION 


The stored-energy function for an ideal per- 
fectly elastic material may be completely defined 
in the following manner. If a cube of the material, 
of unit edge, is subjected to a pure, homogeneous 
deformation in which it is deformed into a 
cuboid of dimensions \; by Az by As, then the 
strain invariants J,, J2, and J; are defined as 


[y=AP+A2?+As3", 
I, = d27A3" +A;"A;" + APA?” 
and 
I3=XPA2A3". (1) 


If the material is isotropic in its undeformed 
state, the energy stored in the deformed cube is 
defined as the stored-energy function for the 
material. It must be a function of J;, Io, and Js, 
provided it depends only on the state of deforma- 
tion of the cube and not on the manner in which 
the deformation is reached. Once the form of this 


2M. Mooney, J. App. Phys. 11, 582 (1940). 
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function is specified, the mechanical properties 
of any body of this material are, in principle, 
completely determinate and are purely a matter 
for calculation, provided such phenomena as 
internal friction and creep are negligible. 

For an incompressible material, the volume of 
the cube considered is unaltered in the deforma- 
tion, so that 


I3=AA27A3" = 1. (2) 
Then W, the stored-energy function, is given by 
W=W(,, I2). (3) 


In general, the form of the function W(J,, Js) 
will be such that when the body is undeformed 
W=0. 


Ill. THE TORSION OF A RIGHT- 
CIRCULAR CYLINDER 


Let us consider a right-circular cylinder of 
radius a and length / of incompressible material, 
for which the stored-energy function W is given 
by (3). It has been shown that, in such a body, 
a pure torsion, in which each circular cross 
section of the cylinder is turned in its own plane 
through an angle ¢(=yz) proportional to its 
distance z from one end, can be maintained by 
the application of surface tractions alone. The 
necessary system of surface tractions is inde- 
terminate to the extent of an arbitrary hydro- 
static pressure, but if the surface tractions 
applied to the curved surface of the cylinder are 
zero, then the surface tractions which must be 
applied to the plane ends have been calcu- 
lated as 


(1) tangential surface tractions 9 acting azimuth- 
ally in the deformed state of the cylinder (these 
effectively act as a torsional couple on the 
cylinder), and 

(4%) normal surface tractions Z. 

8 and Z are measured per unit area of the 
surfaces on which they act and their distribution 
over the plane ends of the cylinder and depend- 
ence on y are given by 


aw aw 
0=24r(—+—-) 
al, aI. 
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and 
ow , ow ow 
Z=— 2a — +29" [ | (—-2—) 
* Ors a ol, dls 
ew ow 
-1|— -+ (3+y?r?)———_ 
ol ,? 01,01» 





ew 
+(2+y’r’) || 4° (4) 
dl,” 


r is the distance of the point considered from the 
axis of the cylinder. OW/dI,, W/dIe, ®W/al,’, 
etc., are, in general, functions of J; and Jz which 
can be calculated if the form of W(J,, Je) is 
known. For the deformation considered, 


I,=1,=3+yr". (5) 


We shall consider here two particular choices 
of the form of the stored-energy function W for 
which Eq. (4) are considerably simplified. 

In the series of papers! mentioned above, par- 
ticular attention was paid to the stored-energy 
function given by 


W=Ci(1i—3), (6) 


when C, is a physical constant characterizing the 
material. Materials for which Eq. (6) are valid 
were described as incompressible, neo-Hookean 
materials. 

The particular interest residing in this choice 
of the stored-energy function rests on the facts 
that 


(1) it is the simplest choice that can be made, 
from the point of view of the formulation of a 
theory of large elastic deformations; 
(iz) it is the form derived by Treloar,* for an 
ideal rubber, on the basis of the kinetic theory 
of large elastic deformations. 

Equation (6) is a particular form, obtained by 
making C.,=0, of the relation 


W=C,\(11—-3)+C2(I2—3), (7) 


which was first suggested by Mooney’? as suitable 
for describing the elastic properties of rubber. 
It is noted that Eq. (7) involves two physical 
constants for the material, C; and C2. 

If the relation (7) is introduced into (4), we 


3L. R. G. Treloar, Trans. Faraday Soc. 39, 241 (1943). 
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a b 


Fic. 1. Distribution of normal surface traction over plane 
ends of cylinder in pure torsion. 


(a) When W=C,(11—3) 
and 
(b) when W =C,(1,-3) + C2(I2—3). 
obtain 

6=2yr(C,+ C2) 
and 


Z= —2Wa’?C.—P* (a — 1’) (Ci —2C2). (8) 


The equations for 8 and Z, when W takes the 
form (6) can, of course, be obtained by putting 
C.=0 in (8). 

We sée that" if the physical properties of the 
material are represented by either Eqs. (6) or (7), 
the force Z which must be exerted at each point 
is proportional to the square of the amount of 
torsion. If Eq. (6) applies, Z must everywhere 
have the nature of a thrust (since C, is essen- 
tially positive), which is distributed over the 
diameter of a plane end as shown in Fig. 1(a). 
If Eq. (7) applies, then Z does not vanish at the 
periphery (where r=a) and provided Cy, is 
positive and sufficiently small compared with C,, 
its distribution will have the form shown in 


Fig. 1(b). 
IV. EXPERIMENTAL PROCEDURE 


The experiments which will be described in 


this paper have as their object the study of the 


distribution of the normal force Z defined in the 
previous section and its dependence on the 
‘amount of torsion. 

In a preliminary qualitative experiment a 
rubber cylinder of pure gum compound, with 
brass plates bonded onto its plane ends, was 
- subjected to a torsional couple by means of a 
lever attached to one of the ends while the other 
end was held in a vise. It was found that unless 
a longitudinal thrust was intentionally exerted 
parallel to the axis of the cylinder, the cylinder 
elongated, its axial cross section before and after 
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Fic. 2. Axial cross section of cylinder before and after 
twisting by a couple. 




















torsion being shown roughly in Fig. 2. This 
simple experiment bears out the mathematical 
conclusion that a pure torsional couple does not 
give rise to a pure torsion. Clearly, if a pure 
torsion is to be obtained, normal compressive 
forces, whose resultant must have the nature of a 
thrust, must be exerted over the plane ends of 
the cylinder. This is in general agreement with 
the theoretical predictions outlined in the pre- 
vious section. 

Consequently, an apparatus was designed in 
which the pure gum cylinder, to whose plane 
ends brass disks are bonded, is restrained from 
elongating. So, when one disk is rotated with 
respect to the other by the application of a 
torsional couple, the longitudinal thrusts re- 
quired for the production of a pure torsion in 
the rubber cylinder are automatically called into 
play. The upper brass disk A contains five 
circular holes, each 3-inch in diameter, varying 
in distance from the center and spaced as shown 
in Fig. 3, in order to have as little mutual inter- 
action as possible. The lower disk B has rigidly 
attached to it a steel disk C. 

The assembly of rubber cylinder and disks A, 
B, and C is held in a robust steel frame D, as 
shown in Fig. 4. The upper disk is rigidly at- 
tached to the frame and the lower disk is capable 
of being rotated in its own plane about its axis, 
by means of a crowbar which is introduced into 
one of a series of holes in the steel disk C. 
Provision is made for holding the cylinder in a 
state of torsion by means of steel pegs, which 
can be introduced into one of the holes E in 
the frame. A scale is also provided for measuring 
the angle through which the lower disk is rotated. 

When this disk is rotated, the rubber bulges 
slightly in the holes in disk A and the amount of 
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this bulging is assumed to be proportional to the 
thrust which would have to be exerted to prevent 
the bulging. 

The bulging represents a small deformation 
superposed on a relatively large torsional de- 
formation. At each of the holes this torsional 
deformation approximates to a state of simple 
shear of an amount depending on the amount of 
the torsion and the position of the hole con- 
sidered. The assumption that the amount of 
bulging is proportional to the normal force tend- 
ing to produce it is equivalent to an assumption 
that the former does not depend greatly on the 
state of simple shear on which the bulging defor- 
mation is superposed. There are certain theo- 
retical reasons for assuming that this is approxi- 
mately valid in the experiments described here. 

In order to produce a simple extension in a 
cylinder superposed on the simple torsion, an 
additional normal surface traction must be 
exerted over each of the plane ends. If the 
stored-energy function for the material of the 
cylinder has the form (6), it has been shown! 
that this additional normal surface traction is 
independent of the amount of torsion and of 
position on the surface of the cylinder. The 
stress distribution in the neighborhood of the hole 
is complex, but for small deformations, it is 
reasonable to assume that the amount of bulge 
in each hole is dependent only on the force 
tending to produce the bulging and not on the 
amount of torsion or position of the hole, except 
insofar as this force is determined by them. This 





Fic. 3. Upper brass disk vulcanized to rubber cylinder. 
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Fic. 4. The torsion apparatus. 


conclusion should be approximately valid if the 
stored-energy function has the form (7), pro- 
vided C2 is small compared with C,. However, 
there is no reason to believe that it should be 
true for every possible form of the stored-energy 
function. 

The amount d of the bulging in each of the 
holes is measured at the center of the hole, for 
various angles of rotation @ of the disk B, by 
means of a dial gauge. The top of the disk B is 
ground flat and used as a reference plane in 
carrying out this measurement. 


V. EXPERIMENTAL RESULTS 


If Eq. (8) is valid for Z, then the relation 
between d and @ should be 


d= K@(2a*C2+ (a? —?r*)(C,—2C2) ], (9) 


where K is some constant of proportionality. 
r is now taken as the distance of the center of 
the hole from the center of the disk. This assump- 
tion is justified if the radius of each hole is small 
compared with that of the cylinder, but it is 
impracticable to make it too small as the values 
of d are thereby reduced. The repeatability with 
which d could be measured in the experiment 
was estimated as +0.002 in. using the same dial 
gauge. The accuracy with which @ could be 
measured in the arrangement employed was +}°. 

The gum compound used had the following 
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Fic. 5. Relation between displacement and the square of 
the angle of torsion. 


composition in parts by weight: 


Rubber 100 
Sulphur 3.0 
Captax 0.8 
Stearic acid 1.0 
Agerite 1.0 
Zinc oxide 5.0 
Tuads 0.08 


The mix was vulcanized for 30 min. at 134°C. 
In order to produce a strong bond between the 


brass plates, a “‘tie gum”’ layer, less than 35” 
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thick, was introduced between the main rubber 


cylinder and each of the brass plates, and a 
cement marketed as ‘“Reanite, No. 3615’ was 
employed. The composition of the ‘tie gum” 
layer was, in parts by weight, 


Rubber 100 
Stearic acid 1.0 
Sulphur 3.5 
Captax 1.0 
Zinc oxide 5.0 
lhermatomic black 40 


The rubber cylinder employed had a radius of 
2.12 in. and length of 1 in. The holes in disk B 
each had a diameter of 3 in. and the distances of 
their centers from the center of the disk are 
given in the following table: 


Hole number 1 2 3 4 5 
r in inches 0 85 1.18 1.47 1.69 


Equation (9) suggests that the results of the 
experiments be plotted as a relation between d 
and @ on linear graph paper. This should yield 
a straight line passing through the origin for 
each hole. 

This is done in Figs. 5 for each of the five 
holes, and it is seen that the linear relation be- 
tween d and & of Eq. (9) is borne out within the 
accuracy of the experiment. 

The slope m of the straight line obtained for 
each hole is the value of K[2a*C2:+(a?—r’) 
X(Ci—2C:2) | for that hole, and if this is plotted 
against (a?—?r’) astraight line should be obtained, 
which passes through the origin if C2.=0 and 
otherwise gives an intercept of 2a°C:K on the 
m-axis. If m is plotted against (a?—r*) as in 
Fig. 6, it is seen that within the accuracy of the 
experiment the relation (9) between d/@ and 
(a*—r*) is again borne out. 

If mo is the intercept of the straight line in 
Fig. 6 on the m-axis and | is its slope, we have 
mo=2a°C2K and b= K(C,—2C2). Whence 
Ci 2b 

=—q?+2. (11) 





Co Mo 


From measurements on Fig. 6, we obtain 


my=10X10~* in. deg.~, 
and 
b=5.65 X10~* in.—! deg. 
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so that 





Mooney? uses the constants C,; and C,2 to 
define three additional constants as follows: 


G= 2(C,+ C2) 
and 


a=I1/G. 


He then shows that the constant G is the modulus 
of rigidity of the material. The constant a@ is 
named the coefficient of asymmetry. From the 
definitions it follows that 


Ci/C2—-1 
a= ete eee 
Ci/Co+1 
Consequently, the experimental results presented 
here lead to a value of 0.75 for the coefficient of 
asymmetry of the pure-gum compound studied. 
This is somewhat larger than the values found 
for two compounds studied by Mooney. Further 
work would be required to establish the variation 
of a with the nature of the compound. 


VI. CONCLUSIONS 


The experiments bear out the general conclu- 
sion of the theory that in order to produce pure 
torsion in a rubber cylinder, normal surface 
tractions must be exerted on the plane ends and 
that these surface tractions vary with the amount 
of torsion and the position on the plane end of 
the cylinder. If the assumption is made that the 
amount of bulging in each hole is proportional to 
the force producing it and does not depend on 
the state of local strain on which the bulging is 
superposed, then the results of the experiments 
are in accordance with the assumption of an 
expression of the form (7) for the stored-energy 
function of the rubber. It has already been 
pointed out by Mooney that a stored-energy 
function of this form describes the linear shearing 
stress vs. amount of shear relation for rubber in 
simple shear, with a fair measure of accuracy, 
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Fic. 6. Relation between slopes of lines in Fig. 5 and 
positions of holes. 


up to reasonably large deformations. Also, 
Treloar* has recently obtained results in agree- 
ment with this for pure, homogeneous strain of 
a sheet of rubber. 

The form (6) for the stored-energy function is 
that derived from the kinetic theory of rubber- 
like elasticity. Load ws. deformation relations 
derived from this have been compared with 
experimental results for simple types of deforma- 
tion and have yielded fairly good agreement, 
depending on the type of deformation employed. 

It appears then that the form (6) for the 
stored-energy function provides a first approxi- 
mation to the stored-energy function for vul- 
canized rubber and that a second approximation 
is provided by the form (7). The differences be- 
tween these forms are more or less accentuated 
depending on the type of deformation studied. 
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Elastic Losses in Some High Polymers as a Function of Frequency and Temperature* 


H. S. Sack, J. Morz,** H. L. Raus, anp R. N. Work 
Cornell University, Ithaca, New York 


The frequency and temperature dependence of the elastic losses of polyvinyl chlorides 
(0° to 80°C, 60 to 1400 c.p.s.) and of natural rubber and butadiene-styrene copolymers (GRS) 
of different relative concentration (— 50° to 35°C, 20 to 3500 c.p.s.) was studied by 3 methods, 
in which the motion of a reed in free or forced vibration is observed. Natural rubber at tem- 
peratures above 20°C, and polyvinyl chloride exhibit frequency independent losses, the mecha- 
nism of which is not yet clear. The results for rubber below 20° and for GRS can be explained 
by a relaxation theory. The behavior of the elastic losses also reflects the existence of a second- 
order transition, below which rotation of chain segments is hampered. 


1. INTRODUCTION 


MONG the theories of elastic losses in solids, 

the so-called relaxation theory has lately 
received considerable attention.' According to 
this theory the reaction of a body to an applied 
strain is not instantaneous but requires a finite 
time which is determined by the molecular con- 
stitution of the material? The time which is 
characteristic of the speed with which molecules 
will move towards their new equilibrium position 
when a stress is set up is called the relaxation 
time. This theory is the analogue to Debye’s*® 


theory of anomalous dispersion of dielectrics, and 


has been discussed for simple and more compli- 
cated models in many publications. It will suffice, 
here, to recall that according to this theory the 


* Presented at the meeting of the Division of High 
Polymer Physics, American Physical Society, January 30 
February 1, 1947. Early results of this investigation were 
presented at the APS meeting in November, 1943. It was 
not possible to publish these results earlier except for a 
short summary (H. S. Sack, J. Motz, and R. N. Work, 
J. App. Phys. 15, 396 (1944).) As they are in direct con- 
nection with the paper presented at the January 1947 
meeting, and in order to avoid duplication in publication, 
the former work, in its essential lines, is included in the 
present contribution. Part of the later work was sponsored 
by the Office of Naval Research. 

More details can be found in Cornell theses of J. Motz, 
R. N. Work and H. L. Raub. 

** Now at the Armour Research Foundation, Chicago, 


Illinois. 


'R. Simha, J. App. Phys. 13, 147 (1942); T. Alfrey and 
P. Doty, J. App. Phys. 16, 700 (1945); W. Kuhn, Zeits. f. 
physik. Chemie B42, 1 (1939). 

? The relaxation theory can also be applied in cases where 
the relaxation is not due to a molecular rearrangement, but 
rather to heat flow. In this case the relaxation time is a 
function of the boundary conditions and therefore of the 
external shape of the body or of the grain structure of the 
substance. In this article we will refer only to relaxation 
phenomena which are governed by the motion of the 
molecules. C. Zener, Metals Tech. [5], 13 (1946). K. 
Bennewitz and H. Rotger, Physik. Zeits. 37, 578 (1936); 
40, 416 (1939). 

3 See e.g. P. Debye and H. Sack; Handb. Rad. (1934); 
Vol.VI, part 1, p. 141. 
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elastic losses—defined as energy loss per cycle 
(often called Q-')—should rise linearly with fre- 
quency w for low frequencies, and decrease as 1/w 
for high frequencies. The shape of the maximum 
in the intermediate region depends on the number 
of relaxation processes and, therefore, on the 
relaxation times that are present. This analogy 
between dielectric and elastic properties has been 
emphasized particularly in a recent publication 
of Kirkwood,‘ in which he applied to elastic 
phenomena a treatment which he formerly had 
applied to dielectric behavior where he calculated 
the distribution of relaxation times on the basis 
of a diffusion process of the single groups of the 
long chain molecules in high polymer substances.® 

So far, little experimental evidence for the ex- 
istence of this relaxation phenomena has been 
found, at least in experiments in which the elastic 
losses under alternating stresses are measured.® 
The aim of the work here presented was to find 
cases in which the relaxation theory would apply. 
High polymer molecules seem particularly suited 
for such a study because these molecules show a 
high degree of mobility, and their motion seems 
to be sufficiently slow so that it can be hoped that 
arelaxation time may be found in a low frequency 
region. It is also known from dielectric measure- 
ments that such relaxation phenomena do exist 
in these substances. 

As a starting-point the series of polyvinyl! chlo- 
rides were chosen because Fuoss’ observed die- 
lectric relaxation in a region around 1000 c.p.s. 
However as will be discussed later, no elastic 


4 J. G. Kirkwood, J. Chem. Phys. 14, 180 (1946). 

5 R. M. Fuoss and J. G. Kirkwood, J. Am. Chem. Soc. 
63, 385 (1941). 

6 The relaxation theory has been tested for nonperiodic 
motions, as for instance in the case of creep, etc. 

7R. M. Fuoss, J. Am. Ch. Soc. 63, 369 (1941). 
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Fic. 1. Diagram of experimental arrangement. (a) and (b) 
coil attachment; (c) arrangement of magnet; (d) optical 
system. 


relaxation was observed. More success was ob- 
tained with elastomers, in particular with syn- 
thetic rubber. The experimental procedure and 
results will be discussed in the following. 


2. EXPERIMENTAL PROCEDURE 


Three different methods for the determination 
of the elastic losses were used : (a) The observation 
of the half-widths of the resonance curve of a reed 
in forced vibration, (b) the measurement of the 
damping properties of a free vibrating body, and 
(c) the measurement of the (elastic) hysteresis 
loop. 

The first and second methods are only appli- 
cable if the damping is relatively low. The last 
method is, in principal, free from this restriction, 
but the interpretation of the measurements is 
easier for high damping. The last method has the 
advantage, that it permits the measurement of 
the elastic losses over a wide range of frequencies, 
without need for changing the sample, while for 
methods (a) and (b) samples of different lengths 
have to be studied, in order to measure the elastic 
losses at different frequencies. 

In all three cases, the sample, in the form of a 
reed approximately 2 mm thick and. 4-15 mm 
wide, and of appropriate length, was clamped on 
one end, while an alternating force could be 
applied to the other end. In the first experiments 
this force was applied by an electrostatic method : 
One side of the reed was covered with a very thin 
metal foil and opposite the reed, separated by a 
few mm, was placed an electrode so that an 
alternating potential could be applied between 
the reed and this electrode. This gives rise to an 
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Fic. 2. Photograph of vibrator-unit. 


alternating force of twice the frequency of the 
applied potential. (If a d.c. potential is super- 
imposed, it is possible to create an alternating 
force of the same frequency.) In latér measure- 
ments a small coil was attached to the end of the 
reed, and the reed was so mounted that this coil 
is in the air-gap of a strong speaker magnet. An 
alternating current set up in the coil then creates 
an alternating force on the reed. 

The motion (vibration) of the end of the reed 
is observed by an optical method. The end of the 
reed protrudes into a light beam. If the reed 
vibrates, this light beam is thus modulated, and 
in falling on a photo-cell creates an a.c. signal 
proportional to the amplitude of vibration. This 
signal is then amplified. Figure 1 shows diagram- 
matically the optical arrangements and the at- 
tachment of the coil to the reed, in (a) for a lateral 
vibration, and in (b) for a longitudinal vibration. 
In the later experiments the whole device was 
made as small and compact as possible in order to 
avoid outside disturbances. Figure 2 reproduces 
a photo of the complete assembly mounted on a 
solid base plate. This picture represents the case 
of longitudinal vibration. For flexual vibrations a 
different clamping block is used. The screws on 
top of the%light-bulb holder serve to adjust the 
light beam so that maximum modulation and, 
therefore, maximum sensitivity are obtained. For 
the study of variations with temperature the 
complete system was put in an automatically 
controlled constant-temperature box. (This ar- 
rangement is, of course, only workable for tem- 
peratures not higher than approximately 70°C.) 

With the compact optical system as shown in 
Fig. 2, it was necessary to define the end of the 
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reed more sharply by attaching to it a thin metal 
edge. The photo-current is first amplified by a 
preamplifier driven by batteries, and then 
by a conventional audiofrequency amplifier. The 
driving potential or current was delivered by an 
audio oscillator, and kept constant over the whole 
frequency range. The frequency could be recorded 
within 1 percent. 

In case (a) (resonance curve), the output of the 
photo-amplifier was read on an a.c. meter giving 
a direct measure of the amplitude of vibration, 
which, when plotted as a function of frequency, 
gives a resonance curve whose half-width is easily 
determined. This method is precise for low or 
moderately high damping. For higher damping, 
the resonance peak becomes very wide and may 
be distorted by higher harmonics or parasitic 
reasonances in the system. 

In case (b) (free vibration), the reed is deflected 
from its rest position by a d.c. potential or cur- 
rent. By opening a switch, the current or po- 
tential is interrupted and the reed returns to its 
equilibrium position executing a damped vibra- 
tion. The output of the photo-amplifier is applied 
to the deflecting plates of an oscilloscope. At the 
time, the switch is opened, starting the vibration 
of the reed, a linear sweep of adjustable speed is 
applied to the other set of the plates of the 
oscilloscope. Thus, on the screen appears the 
picture of the damped vibration, from which the 
logarithmic decrement can be calculated by 
measuring successive amplitudes. Figure 3 re- 
produces two sets of pictures obtained with two 
materials of different constitution. The four 
pictures in each set represent different frequencies 
of vibration as indicated by the time-scale. The 
set at the left represents a material in. which the 
damping increases very strongly with increasing 
frequency, whereas in the set at the right the 
damping is nearly independent of frequency. 
Again, this method is only usable for low 
damping. 

In case (c) (hysteresis loop), the procedure is 
somewhat more complicated. The output of the 
photo-amplifier is again applied to one set of 
plates of the oscilloscope. To the other pair of 
plates a potential is applied which is taken across 
a resistor in series with the driving coil. Since the 
driving force is proportional to the current in the 
coil, this signal is proportional to the force. Thus 
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TYPICAL OSCILLOGRAMS 


ae 





Fic. 3. Examples of damped free vibrations: butadiene- 
styrene, left 50/50, right 70/30. 


there will appear on the screen the elastic 
hysteresis loop of the material (stress-strain 
diagram). The area of this hysteresis loop is a 
measure of the damping qualities of the material. 
In fact, if divided by the square of the amplitude, 
it is proportional to the energy dissipated per 
cycle divided by the total energy of vibration. In 
this method, the two signals must appear on the 
oscilloscope screen without phase differences 
caused by the electric system. As phase differ- 
ences in amplifiers are unavoidable, a phase 
shifter is added to the system. This phase shifter 
is calibrated for each frequency by replacing the 
sample by a metal reed. It is known that for a 
metal like brass the phase shift between the 
vibration and the applied force is near 0 or 180°, 
except for a narrow region around resonance. If 
with this metal reed the stress-strain curve ap- 
pearing on the oscilloscope screen is not a straight 
line, then there exists a phase shift in the elec- 
trical system, and it can be compensated for by 
an adjustment of the phase shifter. The accuracy 
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of this phase adjustment proved to be one of the 
main difficulties of this method, especially for 
measurements at frequencies far away from 
resonance. 

In order to have a relative measure of the 
amplitude, another signal derived from a resist- 
ance in series with the coil is passed through the 
photo-amplifier system. As the current in the coil 
is kept constant, a variation of this signal as it 
appears on the oscilloscope screen is a measure of 
the change in amplification of the photo-amplifier. 

The pattern on the oscilloscope screen is photo- 
graphed, and the area of the hysteresis loop 
measured with a planimeter. This area is then 
divided by the square of the maximum amplitude 
as determined with the help of the calibration 
line mentioned above. 


3. INTERPRETATION OF MEASURED QUANTITIES 


The quantities determined by the three sets of 
experiments are: the half-width of the resonance 
curve Aw/w, (where Aw is the difference in fre- 
quency for two points where the amplitude of 
vibration is 4 the amplitude at resonance fre- 
quency w,) ; the logarithmic decrement 6; and the 
energy W dissipated during one cycle, divided by 
the square of the amplitude y, i.e., W/y?. In 
order to discuss the elastic properties in terms 
of molecular properties, it is necessary, first, to 
relate the measured quantities to material con- 
stants. In the case of forced vibrations this can 
be done by introducing a complex modulus of 
elasticity, e.g., a'complex Young’s modulus for 
longitudinal vibrations, E’+7E”, which might 
be frequency dependent. For small damping 
(Ek /F’<«1) the same complex modulus can be 
used in treating the case of free vibrations. In the 
case of (E’’/E’)<1, the following relations exist 
Aw/w,= V3(E"/E’); 6=2(E"/E’), where E” and 
kK’ have the values at the resonance, or free 
vibration frequency, respectively. W/y* is pro- 
portional to E”’, with a proportionality constant 
which depends, in contrast to the two other 
formulas, on the dimensions of the reed, the 
boundary conditions etc. (For a free-clamped, 
thin reed of length / this constant is r?/8/.) This 
proportionality is only valid for frequencies 
below a certain frequency, which is low for low 
damping and high, for high damping. For fre- 
quencies beyond this limit, deviations from pro- 
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portionality occur that have to be calculated for 
each particular case. The same complication 
arises for the other two methods if the damping 
is no longer small. In all these cases of non- 
proportionality it is necessary, in order to calcu- 
late the higher approximations to know the exact 
mechanism which gives rise to the losses, and 
which determines the frequency dependence of 
the elastic modulus. These calculations were 
made under the assumptions of a thin reed, 
clamped at one end, and acted upon by the ex- 
ternal force on the other end, whose motion was 
observed. The damping mechanism was assumed 
to be of the viscous type, and the frictional force 
was written as K(de/dt), where ¢ is the deforma- 
tion. In this case E’’= Kw. Further calculations 
were made with putting E’’=const., without 
specifying the mechanism that would yield this 
type of loss, since such a mechanism is not yet 
known. The experimental results were then cor- 
rected with the help of the formulas thus derived. 
Since the assumptions are only approximations to 
the real case, no data were used where the 
corrections exceeded 30 percent. 

In the hysteresis loop method, the evaluation 
of W/y* in absolute units is rather cumbersome. 
Therefore the width of the resonance curve was 
measured simultaneously, at least at one temper- 
ature. This value was used to calculate the pro- 
portionality constants. The results were ex- 
pressed directly in terms of Aw/w,, and are 
plotted in these terms in the figures. (E’’(w) is 
then the value taken from the graph multiplied 
by £’ at the resonance frequency and at the 
temperature where resonance is observed, and 
divided by v3.) 


4. EXPERIMENTAL RESULTS AND DISCUSSION 
(a) Generalities 


In all the experiments the amplitude of vibra- 
tion was so small that the maximum strain was 
not bigger than 0.1 percent. Under these condi- 
tions it proved that all the effects are linear. With 
the exciting force varying by a factor 5 to 10 in 
amplitude, no change in the losses was observed. 
Also the vibration remained sinusoidal (no ap- 
pearance of harmonics). Furthermore the hyster- 
esis loop has the form of an ellipse within a very 
good approximation (deviations from the ellipse 
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Fic. 4. Young’s modulus and elastic losses of polyvinyl- 
chloride with 20 percent (A), 30 percent (B) and 40 percent 
(C) plasticizer. 


correspond to an harmonic content of not more 
than 1 or 2 percent). 

The influence of the “tightness” of clamping 
appears to be negligible for the loss measurement. 

The metal foil glued to the reed in the electro- 
static method, and the coil and coil supports used 
in the electrodynamic method, were replaced by 
similar attachments of different weight or geo- 
metric form, and proved to have no influence 
upon the measured elastic losses. 

In most cases measurements were repeated, if 
possible with reeds of different width or length. 
Good reproducibility was observed. The error for 
the losses is approximately +10 percent. 

As will be seen later, very high losses, of an 
equivalent of approximately 100 for Aw/w, were 
observed under certain conditions. Such high 
values have to be interpreted with caution, as 
they mean that all the energy supplied by the 
outside force is dissipated in a very small region 
near the point where the force acts, and it is not 
excluded that the way in which the coil is con- 
nected to the reed has an effect ; also the formula 
used for the calculation will probably receive 


additional correction terms, which cannot be 


evaluated exactly. 


(b) Polyvinyl Chlorides 


Three samples were measured.’ Material A 
contained 80 percent polyvinyl chloride and 20 
percent tricresyl phosphate, B, 70 and 30 percent, 
and C, 60 and 40 percent. All measurements were 


8 We are obliged to Dr. Fuoss, at that time at General 
Electric, for providing us with these samples. 
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Fic. 5. Losses of natural rubber. 


made by the first method (resonance curve). In 
Fig. 4 the results for Young’s modulus and the 
half-widths of the resonance curve are indicated 
as a function of temperature. Young’s modulus 
was calculated from the measured resonance fre- 
quency assuming a free-clamped, thin reed (or a 
clamp-clamped reed) in flexural vibration. By 
choosing different lengths of the reed, and by 
vibrating it as free-clamped or clamp-clamped, a 
series of resonance frequencies are obtained in the 
range of approximately 100 to 1400 c.p.s. The 
values obtained at these different frequencies 
agree within the experimental error with the 
curves drawn in Fig. 4. We can, therefore, con- 
clude that in this frequency range, no dispersion 
occurs. This means that in our experiments we do 
not find a relaxation phenomena as it would have 
been suggested by the dielectric measurements 
made by Fuoss on the same substances in the 
same frequency and temperature range. This 
does of course not necessarily mean that such 
relaxation phenomena are absent, but if they are 
present, they are masked by other processes 
which are of greater importance and yield fre- 
quency independent losses. 

From Fig. 4 it can be seen that Young’s 
modulus decreases with increasing temperature, 
while the losses increase. This can be explained by 
the ‘loosening’ of the structure, and the in- 
creased tendency to flow. It should also be men- 
tioned that static measurements of Young’s 
modulus, at room temperature and with small 
stresses, were made, which yielded from 4 to 10 
times higher values than those obtained by the 
vibration method. This seems to indicate that a 
dispersion zone must be present below the fre- 
quency range here studied. 
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(c) Natural Rubber 


The sample studied was natural gum stock 
with 3 percent sulphur. For temperatures between 
5 and 55°C measurements were made with all 3 
methods, yielding identical resuts ; for lower tem- 
perature, down to —46°C, only the hysteresis- 
loop method proved workable. The frequency 
range was 20 to 900 c.p.s. for the first two 
methods, and 300 to between 1500 and 3500 
(depending on the temperature) for the third 
method. Figure 5 reproduces the results for tem- 
peratures below 4°C, and shows that the losses 
rise linearly with frequency (but the line does 
not go through the zero). Measurements at 
temperatures above 10°C are not included in the 
diagrams. They yield for 20 to 1300 c.p.s., a 
constant value of 0.2 for 22°C, and slightly lower 
for higher temperature. (A slight rise of approxi- 
mately 10 percent would also be compatible with 
the measurements. ) 

Where a comparison can be made, our results 
are in good agreement with those of Stambaugh,’ 
and other earlier investigators. The frequency 
independent losses at higher temperatures can be 
represented, in a formal way, by a constant E”’. 
However, as mentioned, no simple mechanism is 
known that would yield this result. Some specula- 
tions concerning this mechanism were given earli- 
er,!° but no new evidence in favor of or against 
it has been obtained. The losses at lower tempera- 
tures can be represented by the superposition of a 
frequency independent part and a viscous term, 
E"’ = Kw, where K could be called a coefficient of 
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Fic. 6. Temperature dependence of ‘viscosity coefficient”’ 
(slope of loss-frequency curve). 


*R. B. Stambaugh, Ind. Eng. Chem. 34, 1358 (1942). 
10 H.S. Sack, J. Motz, and R. N. Work, J. App. Phys. 15, 
396 (1944). 
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viscosity. The same frequency dependence would 
result from a relaxation phenomenon with meas- 
urements made below the relaxation time (or 
zone of relaxation times). In fact, for low fre- 
quency, the relaxation theory gives a molecular 
picture of viscous flow. Figure 6 shows how this 
coefficient of viscosity increases with decreasing 
temperature. 

Resonance frequencies could be observed above 
—18°C. These values indicate that Young’s 
modulus is rising slightly with decreasing temper- 
ature (20 percent increase from 22°C to —18°C). 


(d) Butadiene-Styrene Polymers 


The samples" were unloaded and contained 1 
percent sulphur, 0.5 percent BLE and 0.3 percent 
tetramethylthiuramonosulfide. The content of 
butadiene to styrene was 70/30; 65/35; 60/40; 
55/45; 50/50. For the first 4 samples, at 22°C, all 
3 methods were used, and gave identical results. 
The last sample at all temperatures, and the 
others below 22° were only investigated by the 
hysteresis-loop method. As an example Fig. 7 
reproduces a set of curves observed for the 60/40 
sample. For higher temperatures than those given 
in the figure, the losses are lower, and the bending 
is less pronounced. In Fig. 6 the slope of the linear 
part (low frequencies) is given as a function of 
temperature for this sample, and for 70/30 and 
50/50. It appears that the temperature variation 
of this slope has the same shape, but is shifted to 
higher temperatures wih higher styrene content. 

The characteristic feature in Fig. 7 as compared 
to the figure for natural rubber (Fig. 5) is the 
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Fic. 7. Losses for 60/40 butadiene-styrene copolymer. 


1! These samples were kindly given to us by the U. S. 
Rubber Company. 
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appearance of a peak in the curve. Again, the 
other samples show a similar behavior, except 
that the temperature at which the bending of the 
curve occurs, in the frequency range here studied, 
is higher for higher styrene content. For the 
70/30 sample it does not yet occur above — 26°C, 
while for 50/50 the bending is already pro- 
nounced at 16°C. This appearance of a peak can 
be interpreted as the result of a relaxation 
phenomenon, though the curves do not corre- 
spond exactly to those calculated with a single 
relaxation time. 

Finally it is of interest to note that while at a 
constant frequency the losses increase first slowly 
and then more and more rapidly with decreasing 
temperature, this trend is reversed at very low 
temperatures (approx. — 40°C for 60/40), and the 
losses decrease again. This fact stands out par- 
ticularly well in looking at the change of the 
resonance curve. In lowering the temperature the 


resonance becomes wider, then is no longer ob- 
servable, reappears again and becomes narrow. 
This characteristic temperature is higher for 
higher styrene content; for 70/30 it could not be 
observed above —50°C, the lower limit of our 
temperature range. This characteristic tempera- 
ture is approximately the second-order transfor- 
mation temperature of these substances.” As it is 
assumed that below the transformation point the 
segments of the chain-molecules ‘‘freeze in,’’ and 
are no longer able to orient as freely as above this 
temperature, we may conclude that the relaxation 
phenomena observed are linked to this segment 
mobility. In agreement with this is the observa- 
tion that the loss-frequency curve below this 
temperature has a different character from that 
above, and does not show the characteristic peak, 
at least not in our frequency range. 

2R. F. Boyer and R. S. Spencer, Advances in Colloid 


Sciences (Interscience Publishers, New York, 1946), Vol. II, 
p. 1 





Electrical Conductivity of GR-S and Natural Rubber Stocks Loaded with 
Shawinigan and R-40 Blacks*: ** 


P. E. Wack, R. L. AntHony, AND E. GutH 
University of Notre Dame, Notre Dame, Indiana 


Measurements have been made of the direct current con- 
ductivity of rubbers loaded with carbon black. Shawinigan 
and Continental R-40 blacks compounded in natural 
rubber and GR-S were studied, and the resistivities were 
determined as functions of time, temperature, concentra- 
tion and elongation. Resistance decreased with time, at 
first very rapidly, then more slowly, approaching an 
equilibrium value. This behavior seems to be independent 
of the type of black used. Temperature coefficients of 
resistance (at 50°C) were positive for Shawinigan stocks, 
negative for samples containing R-40, and tended to 


I, INTRODUCTION 


ARBON black has been used as a reinforcing 
agent in rubber-like materials for many 
-years. In recent years special grades of carbon 
black have been developed which serve the pur- 
* This work was carried out under contract to the Office 
of Naval Research as a part of Contract N6-ori-83, Task 
Order 1. 
** Presented at the American Physical Society Meeting, 
Division of High Polymer Physics, held at Columbia 
University on Jan. 30-31 and Feb. 1, 1947. 
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increase with increasing concentration of black and with 
increasing extension of the sample. At low concentrations, 
R-40 gave higher conductance in GR-S than did Shawini- 
gan. At high loadings, Shawinigan gave the higher con- 
ductance in both GR-S and natural rubber. Resistance 
increased with elongation for all stocks containing R-40. 
The Shawinigan GR-S samples showed an inversion; for 
small elongations the resistance increased on stretching, 
but for higher extensions the resistance decreased. The 
results are interpreted on the assumption that the carbon 
black tends to form chains in the rubber matrix. 


pose of not only reinforcing these materials but 
also of imparting to them some electrical con- 
ductivity. However, our knowledge of the exact 
mechanisms or processes, which are responsible 
for this reinforcing action and greatly increased 
electrical conductivity, is still far from complete. 
Both effects are closely related, and increased 
knowledge of one will aid our understanding of 
the other. 

The reinforcing action of fillers has been dis- 
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cussed in a series of theoretical papers by 
Rehner,! Smallwood,? and Guth.’ A general dis- 
cussion of the effects of various types of carbon 
black in GR-S has recently been given by 
Parkinson.‘ The apparent tendency of the carbon 
black particles to form chains in the rubber, and 
possibly even a network structure at high carbon 
black concentrations has been pointed out by the 
aforementioned authors and others. Parkinson, in 
particular, has noted that, in some instances, 
there is strong evidence that the forces between 
carbon particles are greater than the forces be- 
tween the carbon particles and the rubber. The 
possibility that in some cases the chain of carbon 
black particles may serve as the major reinforcing 
agent rather than the individual carbon black 
particles, was suggested by Guth.* 

Past studies of the electrical conductivity pro- 
duced in rubbers by carbon blacks have been 
confined chiefly to investigations of the effects 
resulting from varying the type of black, the 
particle size of the black, and the concentration 
of the black in the rubber. Cohan and Mackey® 
and Cohan and Steinberg® have shown that the 
specific resistance of a rubber compound at 20% 
volume concentration decreases with increasing 
purity of the black, decreases with increasing 
graphitic structure, and decreases exponentially 
with decreasing particle size. The progressive 
decrease in resistance with decreasing particle 
size had also been noted by Amon and Brown’ 
in a study of the electrical resistance of carbon 
black inks. Hall, Buckley, and Griffith? have 
studied tle dependence of the electrical re- 
sistance upon the concentration and structure 
of the black, using Shawinigan black in natural 


rubber and in GR-S. Bulgin,’ in one of the 


‘J. Rehner, J. App. Phys. 14, 638 (1943); Rubber Chem. 
Tech. 17, 865 (1944). 

7H. M. Smallwood, J. App. Phys. 15, 758 (1944); 
Rubber Chem. Tech. 18, 292 (1945). 

*E. Guth, J. App. Phys. 16, 20 (1945); Rubber Chem. 
Tech. 18, 596 (1945). 

4D. Parkinson, Trans. Inst. Rubber Ind. 21, 7 (1945); 
Rubber Chem. Tech. 19, 100 (1946). 

°L. H. Cohan, and J. F. Mackey, Ind. Eng. Chem. 35, 
806 (1943); Rubber Chem. Tech. 16, 918 (1943). 

® L. H. Cohan, and M. Steinberg, Ind. Eng. Chem. 36, 7 
(1944), 

7F. H. Amon, and O. J. Brown, American Ink Maker, 
November, 1941. 

*R. H. Hall, B. P. Buckley, and T. R. Griffith, Can. 
Chem. and Process Industries 29, 587 (1945). 

*D. Bulgin, Trans. Inst. Rubber Ind. 21, 188 (1945); 
Rubber Chem. Tech. 19, 667 (1946). 
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most comprehensive papers published to date 
on this subject, discussed the dependence of the 
electrical conductivity in conducting rubbers 
upon the concentration of the black in rubber, its 
dependence upon the temperature, and its de- 
pendence upon the extension. Bulgin® employed 
a large variety of blacks and studied the behavior 
of these blacks in GR-S, GR-I, and in natural 
rubber. 

In the work reported here, the conductivity has 
been studied much more completely than in 
any previous work as a function of time, tem- 
perature, extension, and concentration of the 
black in the rubber. The behavior of two types of 
conducting blacks, Shawinigan and Continental 
R-40, in both GR-S and in natural rubber, has 
been studied. The conductivity of the blacks 
themselves was also investigated as a function 
of temperature, under a pressure of approxi- 
mately 2000 Ib. per square inch. 

This more complete study was carried out in 
the expectation that the results might indicate 
the state of aggregation of the black particles in 
the rubber and thus lead to conclusions about 
the relative strength of the carbon black-carbon 
black and the carbon black-rubber interaction. 
Since rubber is an insulator, appreciable con- 
duction in a rubber stock loaded with conducting 
carbon black spheres can take place only if these 
spheres form a chain. Such chains in turn may 
form a lattice or network. Our results can be 
understood on the basis of a simple picture of 
the formation and breaking up of chains of 
carbon black particles. Carbon black particles 
are colloidal in size. To a first approximation, 
we may neglect the forces between them. Then 
we can consider the chain formation as a purely 





TABLE I. 
GR-S stocks Natural rubber stocks 

Parts by Parts by 

weight weight 
GR-S 100.0g Smoked sheet 100.00 g 
Zinc oxide 5.0g¢ Zinc oxide 7.85 g 
Pine tar 5.0g Pine tar — 3.00 ¢g 
Stearic acid 2.0¢ Stearic acid 3.30 ¢g 
Sulfur 2.0g Sulfur 2.81 g 
Captax 15g Captax — 0.743 g 
(Shawinigan or variable (Shawinigan or variable 

R-40) R-40) 


Cured for 70 minutes at 40 
lbs. /in.? steam pressure. 





Cured for 90 minutes at 35 
Ibs. /in.? steam pressure. 
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Fic. 1. Fast stress-strain curves at 30°C for GR-S stocks. 
Shawinigan stocks are denoted by dashed lines and R-40 
stocks by full lines. 


kinetic phenomenon regulated only by chance. 
The probability of formation of a chain becomes 
appreciable only when the loading is sufficiently 
high. Changes in temperature and elongation 
cause a modification of structure of the carbon 
black chains. This will also change with time. 
It should be pointed out that the carbon black 
particles are connected to each other by the 
flexible long rubber molecules. In the unstretched 
state, the rubber molecules are curled up, in 
accordance with the kinetic theory of rubber 
elasticity. The flexibility of the rubber molecules 
will permit almost free translational kinetic 
(Brownian) motion of the carbon black particles, 
up to the limit given by the length of the rubber 
molecules. According to this picture, stretching 
of a loaded stock will limit the region in which 
the carbon black particles can move, and at 
sufficiently large extensions, the positions of the 
carbon black particles are more or less rigidly 
fixed. 


Il. EXPERIMENTAL PROCEDURE 


All stocks were mixed on a small laboratory size 
mixing mill and cured in a small platen press. The 


TABLE II. GR-S stocks. 


Shore 
Volume Tensile hardness 

concentration strength Elongation (30 sec 

Black (%) (Ib /in.?) at break (%) value) 
Shawinigan 40 1270 50 97 
Shawinigan 28 2070 300 74 
Shawinigan 19 1430 750 57 
Shawinigan 15 1080 875 46 
R-40 40 1640 150 97 
R-40 28 2880 675 74 
R-40 19 2820 775 61 


R-40 15 3060 925 50 
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Fic. 2. Fast stress-strain curves at 30°C for natural 
rubber. R-40 stocks are denoted by full lines and Shawinigan 
stocks by dashed lines. 


following base formulae were used for all stocks: 
Four different loadings of Shawinigan and of 
R-40 were employed in the GR-S and natural 
rubber stocks. These were 136.9 g, 80 g, 50 g, and 
36.9 g. These loadings gave volume concentra- 
tions of 40 percent, 28 percent, 19 percent, and 15 
percent, respectively, for the GR-S stocks, and 
39 percent, 27 percent, 19 percent, and 15 per- 
cent, respectively, for the natural rubber stocks, 
the slight difference at the two highest loadings 
being caused by the difference in densities of 
GR-S and natural rubber. The volume concen- 
tration is defined here as the ratio of the true 
volume of the carbon black to the total volume 
of the cured rubber. The carbon black volume 
was calculated from its weight, employing the 
standard density value of 1.80 g/cm*. The 
densities (and volume) of the cured rubber stocks 
were determined by hydrostatic weighing at 
room temperature. It is of interest to note that 
the volume of cured loaded stock is equal to the 
sum of the volumes of the cured gum stock and 
of the carbon black, to within 1 percent or better, 
when the carbon black volume is calculated from 
the density value of 1.80 g/cm*. This was checked 


TABLE III. Natural rubber stocks. 





Volume Tensile - ; 
concentration strength Elongation Shore 
Black (%) (Ib /in.?) at break (%) hardness 
Shawinigan 39 1220 50 99 
Shawinigan 27 1630 150 81 
Shawinigan 19 2060 575 67 
Shawinigan 15 2640 500 66 
R-40 39 1340 150 95 
R-40 27 2060 425 78 
R-40 19 2500 425 70 


R-40 15 2140 550 60 
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Fic. 3. View of oven with door opened to show 
mounting rack. 


by hydrostatic weighing of gum stocks and 
loaded stocks. 

The tensile strength, elongation at break and 
Shore hardness are given in Table II below for 
the eight GR-S stocks, and in Table III for the 
eight natural rubber stocks. 

Figures 1 and 2 show stress-strain curves at 
30°C for the 19 percent and 28 percent concen- 
trations in GR-S, and for the 19 percent and 27 
percent concentrations in natural rubber. It is 
seen that for each loading the R-40 gives higher 
rigidity in both natural rubber and GR-S at 
low extensions than does the Shawinigan. At 
high extensions, the situation is reversed. An 
explanation for the higher rigidity of Shawinigan 
stocks at higher extensions will be suggested in 
connection with the interpretation of resistivity- 
extension curves (cf. Fig. 14). The data shown in 
Figs. 1 and 2 were taken with the fast stress- 
strain machine described in an earlier paper."® 

The rubber was cured in the form of standard 
tensile sheets 6 inches square by ;, inch thick. 
After an aging period of 24 hours or more, the 
samples for investigation were die-cut from these 
sheets. These samples were in the shape of square 
headed dumb-bells, the straight center sections 
having a width of 0.25 inches and a length of 2 
inches. The ends of the samples were placed in 
flat aluminum clamps which covered only the 
enlarged end portions of the samples. The clamps 
served simultaneously as a means of stretching 


10S. L. Dart, R. L. Anthony, and P. E. Wack, Rev. Sci. 
Inst. 17, 106 (1946). 
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Fic. 4. View of oven with door closed. The apparatus 
used for measuring resistivities of the blacks is seen in front 
of the oven. 


the rubber and also as electrical contacts to the 
sample. Resistivities of (stretched) samples were 
measured in the direction of the stretch. The 
coating of the sample ends with Aquadag before 
clamping made no appreciable difference in the 
over-all resistance measurements. After clamp- 
ing, the samples were mounted in an oven by 
means of a rack capable of holding eighteen 
samples, each sample being individually mounted 
in a horizontal Pyrex glass tube. One end of 
each sample was held fixed at one end of the 
glass tube, the other end being attached to a 
copper rod which protruded through the door of 
the oven. The elongation of each sample could be 
adjusted and measured by means of these rods. 
The time required, after a temperature change, 
for equalization of temperature throughout the 
oven, including the interior of the samples them- 
selves, was checked by means of a differential 
thermocouple. This waiting time was approxi- 
mately 15 minutes after reaching the desired tem- 
perature, as read on a mercury thermometer 
centrally located in the oven. Good air circulation 
was maintained by means of a small fan mounted 
in the base of the oven. Resistances were meas- 
ured with the aid of a Leeds and Northrup 
Type H.S. galvanometer and a 2-volt storage 
battery. For low resistance samples a Leeds and 
Northrup Testing Set, Model 5430, was also 
used. Specific resistances were then calculated for 
each sample from the known sample dimensions. 
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Fic. 5. Resistivity versus volume concentration of 
Shawinigan black at 70°C. Open circles denote GR-S 
stocks, and black circles denote natural rubber stocks. 


The measurements on the carbon black were 
made in the oven by inserting the powder in a 
small Bakelite tube and maintaining a pressure of 
approximately 2000 lbs./in.2 on the black by 
means of a piston and tension spring device. This 
pressure, along with considerable tamping down 
of the black while inserting it in the cylinder, 
gave apparent densities ranging from 0.36 g/cm* 
to 0.39 g/cm. 

Figure 3 shows the oven with its door opened, 
and the rack in which the samples are placed. 
Figure 4 is a view of the oven with its door closed 
and shows the copper rods projecting through the 
door. The device employed in measuring the tem- 
perature coefficients of resistivity for the carbon 
blacks is seen in front of the oven. 

Reversal of the direction of current flow during 
the measurements, as well as several checks 
employing 60 cycle current, gave no evidence of 
any important polarization effects, either in the 
rubber or in the carbon black. 


Ill. RESULTS 


A. Dependence of Resistivity Upon Carbon Black 
Concentration and Upon Time 


Initial measurements of the resistivities of 
several samples showed that the resistivity for 
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Fic. 6. Resistivity versus volume concentration of R-40 
black at 70°C. Open circles denote GR-S stocks, and black 
circles denote natural rubber stocks. 


each sample was a function of the time which had 
elapsed after the sample had been cut out and 
mounted in the clamps in an unstretched condi- 
tion. In all cases the resistance decreased with 
time, at first quite rapidly, and then more slowly, 
and appeared to be approaching a stable value. 
This behavior is easily explained using our pic- 
ture of the rubber-black system. As the samples 
were molded, cured, and cut, strains were set up 
in their interior. These strains caused a breakage 
in the carbon black chains, thus leading to an 
increase in the resistance. In the course of time 
kinetic (Brownian) motion of the carbon black 
particles will lead to a re-forming of the chain 
structure and thus to a decrease in the resistivity. 
Any small temporary distortion of the sample, 
for example a small extension followed by an 
immediate retraction to zero extension, produced 
a relatively large increase in the resistance, which 
then decayed in the same manner with time again 
to a semi-equilibrium value. This behavior had 
been noted much earlier by Evans"! and has since 
been pointed out by Bulgin and others. The 
approach to a steady state value of the resistance 
can be greatly accelerated by ‘‘annealing’’ the 
samples at elevated temperatures. This effect 
would seem to indicate that any sudden change in 

R. D. Evans, Paper presented at the meeting of the 


Ohio Academy of Science in 1930 and kindly made avail- 
able to us by W. W. Vogt. 
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strain disrupts existing carbon networks, which 
subsequently tend to reform themselves. In the 
work described here all samples were annealed, in 
an unstretched condition, for a period of 10 hours 
at 70°C after they had been mounted in the oven. 
The value of the resistivity attained at the end of 
this initial annealing period is denoted by po in 
this work. (It should be noted here that this value 
of resistivity is not really constant but rather 
continues to decrease slowly with time. However, 
the important consideration is that the rate of 
change of resistance .with time must be so small 
that it cannot appreciably alter the results of 
some additional experiment.) Preliminary meas- 
urements also indicated wide variations between 
the resistivity values for samples which were 
supposed to be identical. Resistivities for samples 
cut from the same tensile sheets were found to 
vary, in some cases, by as much as a factor of ten. 
This was particularly true for the Shawinigan 
black—GR-S stocks at low carbon concentra- 
tions. The Shawinigan black—natural rubber 
stocks were somewhat more uniform. However, 
the continental R-40 stocks were much superior 
in this respect. At high concentrations the varia- 
tions in resistivity values were much smaller for 
all stocks. 

The dependence of po upon the volume concen- 
tration of the carbon black in the rubber is shown 
in Figs. 5 and 6. Figure 5 shows the effect on 
resistivity of four concentrations of Shawinigan 
in both GR-S and natural rubber, and Fig. 6 
shows the behavior of R-40 in these two rubbers. 
In these two figures each point is an average value 
of a large number of separate determinations on 
individual samples. The vertical line drawn 
through each point shows the range of po values 
for each stock. Comparison of Figs. 5 and 6 shows 
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Fic. 7. Annealing curves at 70°C for 19 percent volume 
concentrations of R-40 (open circles) and Shawinigan 
(black circles) in natural rubber. 
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that at low loadings R-40 gives better conduc- 
tivity (lower resistivity) in GR-S than does 
Shawinigan. At the high loadings, Shawinigan 
gives better conductivity in both GR-S and 
natural rubbers. The R-40 natural rubber stocks 
have higher resistivities than the R-40 GR-S 
stocks. In the case of the Shawinigan black, how- 
ever, the situation is reversed and here the GR-S 
stocks have the higher resistivities. The very 
rapid decrease of resistivity with increasing 
concentration of Shawinigan black at low concen- 
trations, may explain in part the wide variations 
found in the resistivity values for identical stocks 
in this region. Slight variations in concentration 
of the Shawinigan, from sample to sample, might 
thus be expected to produce large changes in 
resistivity in the region of low loadings. Large 
variations in the resistivities of supposedly 
identical rubber stocks loaded with Shawinigan 
black were also found by Hall, Buckley, and 
Griffith.* They suggest that the effect may be due 
to variations in the degree to which the structure 
of the Shawinigan black has been broken down 
during the milling and curing processes. 

To minimize the effects of variations in the 
value of po and thus make possible better 
comparisons between different samples, when 
studying the dependence of resistivity upon ex- 
tension, temperature and time, the relative 
resistivity, p/po, is employed in this work rather 
than the actual resistivity, p. Figure 7 shows the 
decrease in the relative resistivity with time 
during the annealing period for 19 percent volume 
concentrations of R-40 and Shawinigan in natural 
rubber. Figure 8 shows the decrease in this ratio 
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Fic. 9. Decay of resistivity with time at various ex- 
tensions for a 19 percent volume concentration of R-40 in 
GR-S. Temperature = 70°C. 


with time, during the annealing period, for the 
same concentrations of these blacks in GR-S. In 
general, this relative decrease with time was 
found to be larger for the GR-S stocks than for 
the natural rubber stocks, and it became pro- 
gressively smaller for both stocks as the loadings 
were increased. 


B. Dependence of Resistivity of Loaded Stocks 
upon Extension 


At the end of the 10-hour annealing period, the 
samples were stretched to various elongations, 
and the temperature was maintained at 70°C for 
an additional 14 hours. During this period the 
decay of resistivity with time was again observed. 
Figure 9 shows the decrease in the relative re- 
sistivity with time, at various extensions, for a 
GR-S stock containing a 19 percent volume 
.concentration of R-40. These curves have been 
cross-plotted at times of 30 seconds and 14 hours 
to show the dependence of the resistivity upon 
the extension. This dependence is shown in 
Fig. 10. Figure 11 shows the decrease in p/po with 
time, at various extensions, for a 19 percent 
concentration of R-40 in natural rubber. The 
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Fic. 10. Cross-plot of Fig. 9 showing the dependence of 
resistivity upon extension at times of 30 sec (open circles) 
and 14 hours (black circles). 


corresponding dependence of p/po upon the ex- 
tension, at times of 30 seconds and 14 hours, is 
shown in Fig. 12. A very similar behavior was 
observed for the 15 percent and 28 percent 
concentrations of R-40 in GR-S and also for the 
15 percent and 27 percent concentrations of R-40 
in natural rubber. It was not possible to obtain 
data as a function of extension for the 40 percent 
loadings in GR-S nor for the 39 percent loadings 
in natural rubber, because of the extreme stiffness 
of these stocks. 

Figure 13 gives the dependence of p/po upon 
time for various extensions of a GR-S stock con- 
taining a 19 percent volume concentration of 
Shawinigan. A very pronounced inversion in the 
value of the relative resistivity as a function of 
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Fic. 11. Decay of resistivity with time at various ex- 
tensions for a 19 percent volume concentration of R-40 in 
natural rubber. Temperature = 70°C. 
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Fic. 12. Cross-plot of Fig. 11 showing the dependence of 
resistivity upon extension at times of 30 sec (open circles) 
and 14 hours (black circles), 


extension is apparent in the range of low ex- 
tensions. These curves have been cross-plotted at 
times of 30 seconds and 14 hours to show the 
dependence of p/p» upon the extension. This de- 
pendence is shown in Fig. 14. Examination of 
Figs. 13 and 14 shows that for high extensions the 
resistivities are almost immediately reduced to 
values which are very considerably less than the 
values which existed just before the extension. At 
extensions in the vicinity of 15 percent to 
20 percent the resistivities, immediately after 
stretching for the first time, tend to become quite 
large and are extremely sensitive to slight differ- 
ences in extension. Reproducibility of relative 
resistivities in this region was poor. The general 
behavior shown in Figs. 13 and 14 is strongly 
suggestive of a maximum in the resistivity 
versus extension relationship in the region of low 
extensions. This might, perhaps, be caused by a 
general breakdown of the carbon black network 
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Fic. 14. Cross-plot of Fig. 13 showing the dependence of 
resistivity upon extension at times of 30 sec (open circles) 
and 14 hours (crosses). 


structure at these extensions, followed by a 
forced realignment of the carbon chains in the 
direction of the extension at higher elongations 
caused by the orientation of the rubber mole- 
cules. A process of this nature would account for 
the increased conductivity in the direction of the 
extension and might even be accompanied by a 
reduction in the conductivity in the direction 
perpendicular to the stretch. The latter possi- 
bility has not yet been checked, since the equip- 
ment employed in this work was designed to 
measure resistivities (and conductivities) in the 
direction of the extension. 

A probable explanation of the different be- 
havior of Shawinigan and of R-40 at large exten- 
sions may be given. At large extensions of the 
material the rubber molecules most firmly bound 
to the carbon black particles will be subject to 
exceptional high extensions and will make ab- 
normally high contributions to the rigidity. 
Therefore, we would expect carbon black to 
increase the rigidity of the material more when 
it is more firmly bound to the rubber molecules. 
That this is the case with Shawinigan is indicated 
by the decrease in resistivity below the initial 
value at high extensions. However the larger 
size of the Shawinigan particles may also be of 
importance both for the increase of rigidity and 
for the decrease of resistance with increasing 
extension. 

An anisotropic conduction of the type sug- 
gested above has been shown by Habgood and 
Waring” to occur in extruded stocks containing 


2B. J. Habgood, and J. R. S. Waring, Trans. Inst. 
Rubber Ind. 17, 51 (1941); Rubber Chem. Tech 15, 146 
(1942), 
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Fic. 15. Decay of resistivity with time at various exten- 
sions for a 19 percent volume concentration of Shawinigan 
in natural rubber. Temperature = 70°C. 


Shawinigan black. A similar asymetry in graphite 
electrodes was also noted by Buerschaper." 

The behavior of Shawinigan black in GR-S for 
the 28 percent loading was similar to that shown 
in Figs. 13 and 14 for the 19 percent loading. In 
the case of the 15 percent loading of Shawinigan 
black in GR-S, the resistivities of the extended 
samples were too high for measurement with the 
present equipment. A similar inversion in the 
resistivity versus extension relationship was 
found for the 27 percent loading of Shawinigan in 
natural rubber, but was not perceptible for the 15 
percent and 19 percent loadings in natural rubber. 
Figures 15 and 16 show the dependence of relative 
resistivity upon extension and time for the 19 
percent volume concentration of Shawinigan in 
natural rubber. 


C. Dependence of Resistivity of Loaded Stocks 
upon the Temperature 


After the samples had been maintained in 
an extended condition for a time of 14 hours at 
70°C, the temperature was lowered in steps to 
about 30°C and the corresponding values of 
resistivity were noted. The temperature was then 
increased in steps, and readings were taken as 
a function of increasing temperature. In general 
the resistance versus temperature relationship 
-so obtained was not a linear one, and the 
two branches of the curve did not coincide. 
Figures 17, 18, 19, and 20 show the temperature 
dependence of resistivity at zero extension for all 
loadings of Shawinigan and R-40 in both natural 


4 R, A. Buerschaper, J. App. Phys. 15, 452 (1944). 
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rubber and GR-S. In general, the agreement be- 
tween data obtained with decreasing and in- 
creasing’ temperatures becomes better as the 
loadings are increased. This is particularly true 
for the Shawinigan stocks. The temperature 
coefficients of resistivity were determined by con- 
structing tangents to these curves at 50°C, this 


temperature coefficient being defined as 


Ap 

a5o°c =— a 

AT’: pso°c 

The dependence of ase upon carbon black 
concentrations is shown in Fig. 21 for GR-S and 
in Fig. 22 for natural rubber. Each plotted point 
is an average of four separate determinations. It 
is seen that asec is positive for all Shawinigan 
stocks, and increases with increasing loadings. 
With one exception, namely, the 40 percent 
volume concentration in GR-S, the temperature 
coefficient for the R-40 stocks is negative and also 
increases with increasing loadings. Since the 
temperature coefficients at 100 percent concen- 
tration should correspond to those of the carbon 
blacks themselves, the possible trend at higher 
concentrations is indicated in Figs. 21 and 22 by 
dotted line extensions of the curves. 

Figures 23, 24, 25, and 26 show the dependence 
of the resistivity-temperature relation upon ex- 
tension for 19 percent volume concentration of 
both blacks in GR-S and natural rubber. Average 
temperature coefficients over the range from 
30°C to 70°C have been calculated from these 
curves as follows. Let A be the average value of 
p/po at 70°C and B the value of this ratio at 30°C. 
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Fic. 16. Cross-plot of Fig. 15 showing the dependence of 
resistivity upon extension at times of 30 sec (open circles) 
and 14 hours (crosses). 
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Fic. 17. Resistivity versus temperature relationships at 


zero extension for four loadings of R-40 in natural rubber. 


Circles indicate decreasing temperature and crosses in- 
creasing temperature. 
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Fic. 18. Resistivity versus temperature relationships at 
zero extension for four loadings of R-40 in GR-S. Circles 
indicate decreasing temperature and crosses increasing 
temperature. 
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Fic. 19. Resistivity versus temperature relationships at 
zero extension for four loadings of Shawinigan in natural 
rubber. Circles indicate decreasing temperature and crosses 
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Fic. 20. Resistivity versus temperature relationships at 
zero extension for four loadings of Shawinigan in GR-S. 
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Fic. 21. Dependence of temperature coefficient of re- 
sistance upon volume concentration in GR-S. Open circles 
denote Shawinigan and black circles R-40. 
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The dependence of this temperature coefficient 
upon extension for 19 percent volume concentra- 
tions of R-40 and Shawinigan is shown in Fig. 27 
for natural rubber and in Fig. 28 for GR-S. This 
coefficient increases with increasing extension in 
all cases. It is positive for the Shawinigan stocks 
and negative for the R-40 stocks. A similar 
behavior was found for the other loadings. 

These results show that the temperature co- 
efficient of the resistivity of the loaded stocks 
depends strongly upon the concentration and 
upon the extension. The positive coefficient is 
probably connected with the breaking and re- 
forming of carbon black chains. Such a kinetic 
effect will, in general, increase with increasing 
temperature, in contrast to the resistivity of the 
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Fic. 23. Relative resistivity versus temperature at vari- 
ous extensions for a 19 percent concentration of R-40 in 
natural rubber. 


black itself, which will decrease with rising tem- 
perature. The explanation of the negative coeffi- 
cient at lower black concentrations needs further 
consideration. It will also be noted that when a 
run is made first decreasing the temperature 
and then increasing it again to the original value 
the final resistivity is (with one exception: 19 per 
cent concentration Shawinigan in GR-S) larger 
than the starting value. This fact and the pres- 
ence of hysteresis effects in general indicates 
that any disturbances of the material, whether 
thermal or mechanical, tends to break up some 
chains and thus increase the resistivity tem- 
porarily. This effect is simply superimposed on 
the other trends which would be observed if 
measurements could be made on samples in 
complete statistical equilibrium. The decrease of 
the area of hysteresis loops on elongation for 
R-40 and Shawinigan in GR-S may be explained 
as due to the fixation of black particles as ex- 
plained in connection with the effect of elongation 
on resistivity in general. 


D. Dependence of Resistivity of the Carbon 
Blacks upon Temperature 


The temperature dependence of the blacks was 
measured with the aid of the apparatus described 
in Section II, a pressure of 2000 Ibs./in.? being 
maintained on the blacks during the measure- 
ments. The dependence of the resistivity upon 
the temperature is shown in Fig. 29 for both 
blacks. The temperature coefficients of resistivity 
at 50°C, as determined from these curves, are 
—0.0021/°C for R-40 and —0.0006/°C for 
Shawinigan. 
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Fic. 24. Relative resistivity versus temperature at 


various extensions for a 19 percent concentration of R-40 
in GR-S. 


It is interesting to note that even at a pressure 
of 2000 Ibs./in.* the apparent density of the 
carbon black ranged from about 0.36 g/cm* to 
0.39 g/cm’, whereas the apparent density of the 
carbon black in the rubber at the highest loadings 
was about 0.7 g of carbon black per cubic centi- 
meter of cured rubber. It had originally been 
expected that pressures of from 1000 to 2000 
lbs./in.2 would result in a stable value of re- 
sistivity for the carbon black, since a previous 
study of the decrease in resistance of the carbon 
black with increasing pressure seemed to indicate 
a leveling off of resistance values in the region of 
about 1000 Ibs./in.? On the basis of this expecta- 
tion, these resistivity values were at first included 
in the plots of the logarithm of the resistivity 
versus volume concentration shown in Figs. 5 and 
6, as points corresponding to a 100 percent 
volume concentration. This led to an approximate 
linear relationship for the R-40 black in GR-S 
and in natural rubber when log resistivity was 
plotted versus log concentration. However, it has 
since been discovered that the resistivity of the 
black itself depended strongly upon its apparent 
density and that even at constant pressure a wide 
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Fic. 25. Relative resistivity versus temperature at 
various extensions for a 19 percent concentration of 
Shawinigan in natural rubber. 
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Fic. 26. Relative resistivity versus temperature at vari- 
ous extensions for a 19 percent concentration of Shawinigan 


in GR-S. 


variation in apparent densities can be obtained 
by tamping down the black while it is being 
loaded in the cylinder. Since the carbon black 
volume, calculated from the density value of 
1.8 g/cm’, plus the gum rubber volume has been 
shown to be equal to the volume of the cured 
loaded rubber, it follows that the “apparent” 
density of the carbon black corresponding to a 
100 percent volume loading should be 1.8 g/cm‘. 
This density value is much greater than any we 
have been able to produce, and the corresponding 
resistivity value is therefore unknown. This 
strong dependence of resistivity upon apparent 
density of the black has been noted by Benson, 
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Fic. 27. Dependence of temperature coefficient of re- 
sistivity upon extension for 19 percent concentrations of 
R-40 (black circles) and Shawinigan (open circles) in 
natural rubber. 
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Fic. 28. Dependence of temperature coefficient of re- 
sistivity upon extension for 19 percent concentrations of 
R-40 (black circles) and Shawinigan (open circles) in GR-S. 


Gluck, and Kaufmann." They also showed that 
the density of the black measured at a constant 
pressure could be increased by as much as a factor 
of two by various densifying processes, for ex- 
ample by ball-milling. 

In contrast to the large variations in the 
resistivity of the blacks produced by variations in 
pressure and density, the temperature coefficients 
of resistivity were found to be nearly independent 
of pressure and density, provided neither the 
pressure nor the density was too low. Accordingly 
these values have been included in Figs. 21 and 22. 


IV. SUMMARY 


The behavior of two types of conducting 
blacks, R-40 and Shawinigan, has been studied in 
natural rubber and in GR-S. A series of four 
loadings of each black in each type of rubber was 
employed. These loadings ranged from about 40 
grams of black per 100 grams of raw rubber to 


4G. Benson, J. Gluck, and C. Kaufmann. In paper 
presented before the Toronto Meeting of the Electro- 
chemical Society, October 16 to 19, 1946. 
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140 grams of black per 100 grams of raw rubber, 
the corresponding true volume concentrations 
ranging from 15 percent to 40 percent. The 
resistivities of these stocks were studied as 
functions of time, temperature, extension, and 
loading. The following results were observed: 


A 


Fast stress-strain curves show that in both 
natural rubber and GR-S, the R-40 black gives 
the better reinforcement at low extensions, but 
that Shawinigan gives the better reinforcement 
at the high extensions. 


B 


Any small strain produces an increase in re- 
sistivity which then decays with time in a nearly 
exponential fashion. In particular the process of 
cutting out the samples and mounting them in 
clamps produces such strains. Annealing the 
samples at elevated temperatures accelerates the 
approach to an initial steady-state resistance 
value. Accordingly, in this work all samples were 
annealed for 10 hours at 70°C in an unstretched 
condition. The resistivity value attained at the 
end of this time is denoted by po. 


C 


At the lowest concentrations, R-40 gave better 
conductivity in GR-S than did Shawinigan, the 
log resistivity versus concentration curve being 
very steep for the Shawinigan black in this region. 
At the high loadings Shawinigan gave higher 
conductance in both GR-S and natural rubber. 
The R-40 GR-S stocks were better conductors 
than the R-40 natural rubber stocks. For 
Shawinigan black, however, the natural rubber 
stocks were better conductors than the GR-S 
stocks. 


D 


For the R-40 stocks, stretching a sample pro- 
duced an increase in the resistivity, which then 
decayed with time while the sample was held at 
a fixed extension. The Shawinigan stocks behaved 
differently, however. In the case of the 19 percent 
and 28 percent concentrations in GR-S and also 
the 27 percent concentration in natural rubber, an 
inversion was found in the resistivity versus ex- 
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tension relationship. The resistivity increased 
with increasing extensions up to about 20 percent 
extension and then decreased for higher ex- 
tensions. For sufficiently high extensions, the 
resistivities were immediately reduced to values 
which were very considerably less than the 
unstretched values. All resistivities were meas- 
ured along the length of the sample and hence in 
the direction of the extension. 


E 


The dependence of resistivity upon tempera- 
ture was observed for all stocks over the tem- 
perature range from 30°C to 70°C. In general, the 
average temperature coefficient, calculated at 
50°C, was positive for the Shawinigan stocks and 
negative for the R-40 stocks. It increased with 
increasing carbon black concentration and also 
increased with increasing extension of the sample. 


F 


Temperature coefficients of resistivity were 
also measured for the two blacks at 50°C. The 
blacks were compressed under a pressure of about 
2000 Ibs./in.2 This maintained 
during the measurements, the ‘“‘apparent”’ densi- 
ties ranging from 0.36 to 0.39 g/cm*. The ob- 
served temperature coefficients were —0.0021/°C 
for the R-40 black, and —0.0006/°C for 
Shawinigan black. 


pressure Was 


‘ 


Our results can, in general, be understood on 
the basis of a simple picture of the formation 
and the breaking up of chains of carbon black 
particles. Thermal and mechanical disturbances 
both lead to a breakage of the chains and thus 
to an increase in the resistance. In the course 
of time Brownian motion of the carbon black 
particles will re-form the chain structure and 
thus decrease the resistivity as a function of time. 
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Fic. 29. Dependence of carbon black resistivity upon 
temperature at a pressure of 2000 lbs. /in.? 


This explains the decrease of resistance with time 
of freshly cured and cut and of stretched samples, 
the increase of resistivity on stretching in general, 
the fact that the ascending branch of the resis- 
tivity-temperature curve lies above the descend- 
ing branch and a positive temperature coefficient 
of resistivity. In the case of Shawinigan the 
greater increase of rigidity and the decrease of 
resistivity (below the initial value), both at high 
extensions, is probably determined by the greater 
strength of the rubber-carbon black bond and 
the larger size of the black particles as compared 
with R-40. 

It is a pleasure to thank Professor Hubert M. 
James for frequent consultation and many valu- 
able suggestions in connection with this study. 
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Elastic Properties of Cork 


II. Stress-Temperature Relationship of Compressed Cork 


S. L. Dart AND EUGENE GUTH 
Polymer Physics Laboratory, University of Notre Dame, Notre Dame, Indiana 


(Received October 14, 1947) 


In the first paper of this series it was shown that for cork, held at constant compression, the 
stress as a function of compression, temperature, and time was separable into two factors, one 
factor representing the load-compression curve and the other describing the decay of stress with 
time. In the present paper the dependence of the load-compression curve upon temperature is 
studied. After the sample has relaxed at constant temperature and compression for a time 
sufficiently long so that the change of stress with time may be neglected, stress-temperature 
curves were obtained. These curves showed a linear dependence of the stress upon the tempera- 
ture. From the stress-temperature curves the resolution of the stress into its additive com- 
ponents due to internal energy and entropy was carried out. The physical significance of this 
resolution is discussed. The experiments were carried out with a refined version of the apparatus 
described in the first paper. In particular, the apparatus was corrected for rate of thermal 


expansion. 





INTRODUCTION 


N a previous paper! an account was given of 

the stress relaxation of cork as a function of 
compression, temperature, and time. In particu- 
lar it was shown that the stress decayed to zero 
in about 1000 years at room temperature and in 
about 1000 minutes at 200°C, independent of 
compression. This independence of the decay 
time on compression was also checked at 80°C 
and at 130°C. It was also demonstrated that the 
compressive stress Z as a function of compression, 
temperature, and time was separable into two 
factors, 


Z(e, T, t)= Fle, T)-G(t, T), (1) 


where the factor F(e, 7) is essentially the load- 
compression curve of cork and the factor G(t, T) 
describes the decay of stress with time and 
temperature. In the previous paper it was proved 
that G(t, T) has the form 
G(t, T) =a—b logt (2) 
where a and 6 depend on T only. In particular 
a and b were thrown into the form 
a=1 
b=8/(a—T) 
giving as a final equation 
Z(e, T, t)=A(e, T)[1—[B/(a—T) ] logt]. (3) 
1 Dart and Guth, J. App. Phys. 17, 314-318 (1946). 
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However, the dependence of A(e, T)=F(e, T) 
upon JT was left open. The present paper fills 
in the gap and brings the fundamental stress- 
temperature relationship for cork. 


EXPERIMENTAL METHOD 


The apparatus used was described in the 
previous paper. It consisted of a lever arm for 
compressing the sample which was automatically 
kept balanced through a chainomatic drive. 
Through this drive the force was recorded as a 
function of time on a clock-driven drum. 

The apparatus was carefully rechecked for 
possible errors in the thermal expansion correc- 
tion and was also corrected for rate of thermal 
expansion. This latter is very important in 
connection with the measurement of stress as 
a function of temperature at constant com- 
pression, since the compression must be main- 
tained constant during the temperature change 
as well as during the temperature equilibrium. 
The procedure used was patterned after that of 
Anthony, Caston, and Guth.? That is, the sample 
was allowed to relax! until the rate of change of 
stress was small enough to make the stress- 
relaxation negligible during the subsequent stress- 
temperature determination. The preliminary re- 
laxation time used was 1000 minutes. After this 
~ 2R. L. Anthony, R. H. Caston, and E. Guth, J. Phys. 


Chem. 46, 826 (1942). Reprinted in Rubber Chem. Tech. 
16, 297 (1943). 
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preliminary relaxation period the temperature 
was lowered in steps and the corresponding stress 
changes were noted. Using this procedure one 
can find the temperature behavior of the stress 
at temperatures less than the relaxation tem- 
perature. 

Figure 1 shows the results of such an investiga- 
tion at the compressions indicated and with the 
common relaxation temperature of 80°C. It is 
seen that within the experimental error the stress 
depends linearly on the temperature. Thus, we 
may express this by the equation 


Z(e, T)=M(e)+N(e)T, (4) 


where M(e) and N(e) are functions of compression 
alone. We will defer the discussion of the sig- 
nificance of this relation to a later section. 

This then gives us the stress as a function of 
compression and temperature, exclusive of time 
effects. That is, this Eq. (4) gives us the tem- 
perature dependence of the factor A(e, 7) or 
F(e, T) in Eq. (3). Before we write down an 
equation expressing all the experimental results 
presented, we must recognize that we have 
measured time effects and temperature effects 
separately, i.e., time effects were measured at 
constant temperature and vice versa. Let us call 
the constant relaxation temperature 7* and the 
temperature which is changed, 7. 7™* is really a 
variable separate from JT because in order to 
specify the state of the system one must specify 
not only the actual temperature 7, but also the 
temperature J* at which the relaxation took 
place. This situation was discussed in more 
detail by Guth, Anthony, and Wack* for the 
case of rubber in extension. They showed, experi- 
mentally, that certain ratios (cf. Eqs. (13a, 
b, c, and d) of their paper) are constant. The 
constancy of these ratios implies the validity of 
‘generalized similarity and separability’ ex- 
pressed by the equation 


Z(e, T, T*, t) = F(e, T)-G(T™*, t). (5) 


It seems reasonable to assume that “‘ generalized 
similarity”’ will also hold for cork, though this 
was not checked directly. Taking into account 
(3) we can then identify more explicitly the 


3E. Guth, P. E. Wack, and R. L. Anthony, J. App. 
Phys. 17, 347 (1946). 
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Fic. 1. Stress-temperature curves for cork previously 
relaxed at T*=80°C for 1000 minutes. Circles indicate 
points obtained on lowering the temperature. Crosses are 
the points obtained on raising the temperature. 


factors F and G. 


F(e, T)=M(e)+N(e)-T; 
G(T*, t)={1—[8/(a—T*) ]logt}. (6) 


Then our generalized equation of state takes the 
form 


Z(e, T, T*, t)=[M(2)+N(OT] 
x }1—[B/(a—T*) ] logt} (7) 
= F(e, T)-G(T*, t). 


DISCUSSION OF EXPERIMENTAL RESULTS 
A. Stress Relaxation 


The discussion so far has been concerned solely 
in presenting the experimental results of our 
investigation of the elastic properties of cork. 
We shall now try to give a more fundamental 
discussion of these results in the light of present 
day theories. Consider first, the stress relaxation. 
The general equation for the time dependency 
of stress as given by Boltzmann has been general- 
ized for the non-Hookean case, by Guth, Wack, 
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and Anthony.’ This may be written as follows: 


t 
Z(t) = Fle)]—8f o(t—t')Fle(t’) jdt’, (8) 
where #’ is a small constant, F= F(e) gives the 
stress-strain relationship, and 6’¢(t—?t’) is the 
memory function of Boltzmann. For relaxation 
we have e(¢) = e(t’) =0 for <0 and e(t) =constant 
for >0. Thus, Eq. (8) becomes 


Z(t) = F(e)L1—B'Y(d) J, (9) 


where 
t 
vin= f o—rar. (9a) 
0 


It is to be noticed that Eq. (9) is identical in 
form with the experimental Eq. (3). Thus, we 
can describe the time dependency by means of 
the generalized Boltzmann theory. Also, for 
the case of cork, comparison of (3) and (9) 
shows that 


A’y(t) =(8/(a—T) | logt, (10) 
while (9a) shows that 
o(t—t’) =1/(t-?t’). (10a) 


A relation of the type (10) holds for many ma- 
terials in a limited range of ¢ values. The re- 
markable result for cork is the validity of the 
relation (10) over a very wide range of ¢ values. 
As a matter of fact, (10) is valid for cork in 
the range ¢=1 min. We plan to discuss the 
theoretical significance of this fact at a later 
time. 

If the stress-time curve is an exponential in f, 
we have 


Z=Z,exp(—K’'t), (11) 
where K’ is given by 
K’ = (kT /h)e~ 48/87, 


Here AH is the activation energy and k7/h is 
the rate factor as given by Eyring.‘ This equation 


is a straight line on a logZ vs. ¢t plot. Such a. 


plot was tried on a relaxation run at 200°C and 
4.86 percent compression and the result was far 
~ from linear. However, activation energies were 
calculated using a straight line portion of the 
curve and also using the time for the stress to 
~~ 4Cf. for instance, S. Glasstone, K. Laidler, and H. 


Eyring, The Theory of Rate Processes (McGraw-Hill, New 
York, 1941). 
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drop to a value 1/e times the first value. These 
activation energies were 38.0 and 36.8 kilo- 
calories per mole, respectively. These results are 
slightly larger than the corresponding values for 
rubber-like materials, which have activation 
energies of about 32 kilocalories per mole. 

In Eq. (3), Fle, 7) designates the value of the 
stress at ‘=1 min. Expressing this fact more 
explicitly, Eq. (7) may be written 
Z(e, T, T*, t) = Fe, T) 

X {1-—[8/(a—T*) ] logt/to}, (12) 


where f9=1 min. For a<7* Eq. (12) was not 
established experimentally and it cannot be 
extrapolated to this range without further con- 
siderations. At any rate theoretical arguments 
show that a logarithmical relationship of the 
type of Eq. (12) will not hold for very small 
times, i.e., for times much smaller than 1 min. 

In the first installment of this series' the time 7 
needed for the stress to decay to zero, for a 
sample held at constant compression, was shown 
to be independent of the compression strain. To 
show the dependence of 7 upon the relaxation 
temperature 7*, a plot of 7* versus logr was 
used. The curve so obtained was a straight line. 

Instead of a 7* vs. logr plot, a logr vs. 1/T* 
plot may also be used. A straight line does not 
fit in the latter representation as well as the 
former one, but still gives a reasonable fit, if 
one considers the unavoidable experimental 
errors. We then have the relation 

1/T* =a+b logr. 
Solving for 7 we obtain 
t=rT, exp(V/RT*), (13) 
i.e., an infinite decay time at absolute zero, which 
is reasonable. The constants in Eq. (13) have 
the values 
To = 1.6-10-* min., 


"=9 kilocalories per mole. 


B. Stress-Temperature 


Concerning the stress-temperature relation- 
ship, we can proceed as follows: 

We start from the well-known energy equation 
of thermodynamics 





aU oP 
( ) -1(—) _P. (14) 
OV T oT V 
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Using the Maxwell relation 


oP as 
(2)-($), 
oT V OV T 

we can write 

aU as 

( ) -7( ) —P. (16) 

OV/ + OV/ + 


Applying this to our case of unilateral com- 
pression we must substitute 1—e for V, and Z 
for P, since we plot actually the pressure or the 
negative of the tension stress. This gives 


aU as 
-(3),--€2) 2 
de T de T 
aU as 
z~(—) -1( ) =Zy+Zs, (17) 
de T Oe T 


where Zy is the contribution of the internal 
energy to the stress, Zs is the entropy contribu- 
tion, and Z is, of course, the total stress. (Actually 
the experiments are carried out at constant 
pressure and constant compression. Nevertheless, 
Eq. (17) is valid in sufficient approximation, 
when derivatives are taken at constant tempera- 
ture and constant pressure.) Equation (17) has 
the same form as our experimental Eq. (4). 
Thus, using the curves shown in Fig. 1 we can 
calculate Zy and Zs as a function of e. In doing 
this one must remember that the T in Eq. (17) 
is absolute temperature. The results of such a 
calculation are shown in Fig. 2. Here one can 
see that the entropy is responsible for most of 
the stress in the first two regions, i.e., for the 
stiffness of the sample at small compressions and 
during the breakdown of the structure. However, 
when the actual material of the cell walls begins 
to be compressed, at 60 percent or so, the internal 
energy contribution becomes appreciable. 

The negative internal energy contribution in 
the first range of (small) compressions is at least 
partly caused by the thermal expansion of the 
uncompressed cork. This is similar to the be- 
havior of rubber at small extensions. The internal 
energy contribution is partly caused by the 
thermal expansion of the unstretched rubber. 
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Fic. 2. Resolution of the stress of Fig. 1 at T=25°C 
into its components Zy and Zs due to internal energy and 
entropy, respectively. 


The large contribution to the stress due to 
entropy is presumably connected with the bend- 
ing of the cell walls. The material of the cell 
walls is likely to consist of long chain molecules 
with at least some degree of free rotation. On 
this basis a large entropy contribution can be 
understood. The increase of internal energy, on 
further compression, is presumably caused by 
intermolecular forces coming into prominence as 
compression increases. 
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Elastic Properties of Cork 


III. Hydrostatic and Ordinary Load-Compression Curves for Cork 


S. L. Dart,* Howarp A. Rosinson,** AND EUGENE GuUTH* 
(Received October 14, 1946) 


Load-compression curves were taken for natural cork. The S-shape of the curves is discussed 
and a theoretical fit is made for the high end of the curve, assuming an initial stiffness and an 
incompressible part of the cork material. A comparison is made between repeated compressions, 


with and without a recovery period of time. 


Pressure-volume curves were also taken and their significance discussed in connection with 
the load-compression curves leading to an evaluation of Poisson's ratio. 


ETERMINATION of the load-compression 
curves of cork under various experimental 
conditions is of great theoretical and practical 
interest. The load-compression curves for cork 


G 


























oY 


Fic. 1. Apparatus for hydrostatic compression of cork. 
T: high pressure tank; G: gauge and reducing valve; 
S: scale; C: glass capillary tube; A : compression chamber. 











* Polymer Physics Laboratory, University of Notre 
Dame, Notre Dame, Indiana. 

** Research Laboratories, Armstrong Cork Company, 
Lancaster, Pennsylvania. 
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show a peculiar S-shape. The proper interpreta- 
tion of this S-shape should throw considerable 
light on the structure of cork. A knowledge of the 
structure of cork, in turn, will facilitate the design 
of corky materials for practical applications, and 
enable one to predict their behavior under a 
great variety of conditions. 

In this paper an attempt is made to interpret 
the S-shaped load-compression curve of cork as 
the result of the cellular structure and of the 
two-phase structure of cork, namely, the corky 
material and the air trapped in the cork cells. 
Load-compression curves were obtained in two 
ways: in the ordinary way, and second, by 
measuring the change in volume caused by 
hydrostatic pressure. 

These two methods lead actually to somewhat 
different load-compression curves; their mutual 
relation will be explained in this paper in the 
section on hydrostatic curves. 


DESCRIPTION OF THE EXPERIMENTAL METHOD 
A. Ordinary Load-Compression Curves 


These curves were obtained by means of the 
apparatus described in a previous paper' (cf. 
Fig. 1 there). The same sample was used to 
obtain the complete load-compression curves up 
to 70 percent compression by gradually increasing 
the load, and reading the corresponding com- 
pressions on the dial gauge. 


B. Hydrostatic Volume-Pressure Curves 


These curves were obtained on an experi- 
mental arrangement shown in Fig. 1. A cork 


! Dart and Guth, J. App. Phys. 17, 314 (1946). 
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Fic. 2. Load compression and extension characteristics of 
cork and of GR-S tread. The crosses indicate the curve 
obtained by repeating the loading a second and a third 
time on the same sample. 


sample of about six cubic centimeters volume is 
placed in the compression chamber which is then 
bolted shut. Water is then forced into the 
stopcock at the bottom until it rises near the top 
of the glass tube and scale. The stopcock is then 
closed and any air bubbles are made to rise to the 
top of the tube by tapping. Care must be taken 
to get all air bubbles out as they would cause an 
error in the volume-pressure curve. A zero 
reading on the scale is taken of the water level, 
and the whole unit is connected through a 
reducing valve to a high pressure tank (of nitro- 
gen in our case). The pressure is then increased in 
units of 5 p.s.i. every 30 seconds, and the corre- 
sponding water-level readings on the scale are 
taken at the end of each 30-second period. This 
procedure was used in order to simulate as 
closely as possible the procedure used in taking an 
ordinary load-compression curve. The pressure 
was read from the reducing gauge and the volume 
change was calculated from the scale readings 
(knowing the diameter of the glass tube). It was 
interesting to note that there was a great deal of 
flow of the cork during a run, as is the case during 
an ordinary load-compression run. This implies 
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that the flow is caused by the cork structure 
itself and not by the escape of air from the cells, 
which would be impossible in the hydrostatic 
case. 


ORDINARY LOAD-COMPRESSION CURVES 
FOR CORK 


Figure 2 represents load-compression curves 
obtained at the Armstrong Research Labora- 
tories. Up to the breaking point the extension 
branch of this curve is also given for cork. The 
crosses indicate the load-compression curve for 
cork obtained by repeating the loading a second 
and a third time on the same sample. Twenty-four 
hours elapsed between consecutive loadings. 

Figure 3 shows a somewhat more extensive 
study of the repeated load-compression curves 
carried out at Notre Dame. The consecutive 
curves were taken one after the other without 
giving time for the sample to recover. Complete 
hysteresis loops were obtained. 

Figures 2 and 3 reveal a series of interesting 
features important for the interpretation of the 
load-compression curves. 

The extension branch is monotonic. The tensile 
strength and breaking elongation are low. They 
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Fic. 3. Load-compression curves for repeated compressions 
including hysteresis loops. 
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Fic. 4. Load-compression curve of cork. Crosses indicate 
theoretical Eq. (1) of text. 


are chiefly determined by the cuts or cracks 
invariably present in cork, and the propagation 
or growth of these cracks. The tensile strength 
may be compared with the tear resistance of a 
compact material like rubber, the tear being 
natural and unavoidable for cork. 

The compression branch shows a characteristic 
S-shape. On repeated compression the S-shape 
tends to be washed out and the curve approaches 
somewhat the monotonic behavior shown by 
GR-S tread.2 This is shown with particular 
clarity in Fig. 3, but to a lesser degree in Fig. 2 
because of the recovery of the sample between 
consecutive loadings. 

A simple explanation of the S-shape may be 
based on the known cellular and two-phase struc- 
ture of cork: the corky material (including various 
organic substances) constituting the cell walls 
and the air trapped in the cells. We notice that the 
knee of the compression curve occurs roughly at 
a compression corresponding (on inversion at the 
origin) to the breaking elongation, and for a load 
somewhat smaller than the negative of the tensile 
strength. On repeated stretching the knee occurs 





2 Synthetic rubber compound containing carbon black as 
used for the tread surface of tires. 
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for a lower critical load. The equal order of 
absolute magnitude of the critical load and com- 
pression of the knee and of the tensile strength 
suggest that a certain breakdown of the cells 
may have taken place at the knee. The difference 
in the lowering of the knee as shown in Fig. 2 
with recovery period between repeated loadings, 


and in Fig. 3 without a recovery period, indicates 
that the breakdown of the cells is partially 
reversible. 

For the part above the knee of the compression 
curve we have tried a preliminary quantitative 
representation. That part of the compressive 
stress Z may be represented fairly accurately by 
the equation 


(Z —a)(L—lI») =const; (1) 


L is the relative length, @ and J) are constants. 
The quantity a should be independent of the 
compression, and represents the critical stress 
necessary for the partial destruction of the cork 
structure (not the cells). 

The plateau after the knee is interpreted as 
largely caused by the failure of the resistance of 
the cork (through buckling of the cell walls for 
instance). The upward turn is interpreted as a 
compression of the almost incompressible part of 
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Fic. 5. Load-compression loop for balsa wood. 
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cork. Of course, the van der Waals 6 of the air 
comes in also. The constant on the right side of 
Eq. (1) should be of the order RT. According to 
this equation, Z should be a linear function of the 
temperature at constant compression (since the 
constant on the right side of Eq. (1) should be 
proportional to 7°) 


Z=M(L)+N(L)-T. (2) 


This conclusion was checked in the second paper 
of this series by our relaxation method. Equation 
(1) checks well with a load-compression curve, as 
shown in Fig. 4. However, the true significance of 
this apparent agreement has not been established. 
For some samples the plateau above the knee 
was practically straight. The possible significance 
of this finding, at some variance with Eq. (1), 
needs further study. 

It is interesting to note that balsa wood, when 
compressed perpendicular to the grain, also 
shows a pronounced S-shaped load-compression 
curve, with an almost horizontal plateau, as 
shown in Fig. 5. The large set (46 percent) after 
compression disappears almost completely when 
the sample is soaked in water. 


LOAD COMPRESSION CURVE OF HIGHLY 
COMPRESSED CORK 


In connection with the interpretation of the 
S-shaped load-compression curve it was inter- 
esting to see how a high hydrostatic pressure 
would change the shape of the load-compression 
curve. At our request Dr. P. W. Bridgman* was 
kind enough to apply to cork a pressure of 29,000 
kg/cm? in isopentane. The initial dimensions of 
the sample on April 11, 1945, were 0.324” 0.319” 
0.312”. After exposure to high pressure the 
dimensions were 0.297" X0.301” X0.283”. On 
lebruary 6, 1946, the dimensions were 0.302” 
x 0.306” X 0.287". In view of the magnitude of 
the hydrostatic pressure employed in the test, it 
is remarkable that the permanent set of the 
linear dimensions was only about 8 percent on 
the average immediately after the exposure. 
After ten months’ rest the permanent set was 
somewhat less. 

In Fig. 6 the dotted line shows the load- 
compression curve taken in June, 1946, on the 


* Jefferson Physical Laboratory, Harvard University, 
Cambridge, Massachusetts. 
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Fic. 6. Comparison of load-compression curves for cork; 
dashed line shows curve for sample hydrostatically com- 
pressed by P. W. Bridgman, applying a pressure of 29,000 
kg/cm?. Full line shows behavior of a sample which was not 
previously compressed. 


same piece of cork which had been subjected to 
the high pressure by Dr. Bridgman on April 11, 
1945. The full line shows the load-compression 
curve of Fig. 2, which was obtained on a piece of 
cork not previously compressed. It is seen that 
the initial part of the two curves is almost the 
same. We interpret this initial part as resulting 
from the initial bending of the cell walls. The 
comparison of the two curves exhibits nicely the 
fact that the cells were not crushed by the high 
hydrostatic pressure. The plateau of the curve 
for the pre-compressed sample is, however, much 
less developed than in the case of the previously 
uncompressed sample. We attribute this plateau 
to the failure of the resistance of the cork cells. 
The comparison of the two curves shows that the 
hydrostatic compression must have driven a 
large amount of the air out of the cells. The steep 
part of the two curves has about the same slope, 
as one would expect. 


HYDROSTATIC CURVES 
It is clear from the description of this method 


in the first section that hydrostatic compression 
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Fic. 7. Comparison of hydrostatic (full line) and ordinary 
load-compression (dashed line) curves for cork. 


leads to a pressure-volume relationship, different 
from the ordinary load-compression relation. The 
connection between these two types of relations 
may be exhibited simply and quantitatively, 
referring to the elements of the theory of elasticity. 
The slope at the origin of a load-compression 
curve Z=Z(e) is Young’s modulus 


dZ/de=E. (3) 


Similarly, the slope at the origin of the pres- 
sure-volume relation, P=P(V), is the bulk 
modulus K 


VdP/dV=K. (4) 
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Now according to the theory of elasticity, 
Young’s modulus, E, and the bulk modulus K are 
connected as follows: 


E=3K(1—2y), (S) 


where yu is Poisson's ratio (ratio between lateral 
contraction and longitudinal extension). Knowing 
E and K Eq. (5) enables us to obtain uy. 


u=43(1—E/3K). (6) 


Figure 7 shows a hydrostatic pressure-volume 
change curve with the corresponding load-com- 
pression curve for cork. When inserted into 
Eq. (6), comparison of the slopes of the two 
curves yields 
pw=0.04. (7) 
Another method of obtaining yu is to measure 
the change of volume on ordinary compression. 
According to the theory of elasticity the following 
relation holds: 
AV/V=e(1—2y). (8) 


This relation permits one to compute u if AV/V 
is measured for a given compression ¢. This 
method, however, was not used in the present 
work. 
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A Cold Cathode Rectifier 


WILLARD H. BENNETT* 
(Received October 21, 1946) 


A cold cathode rectifier is described in which corona discharge is used at atmospheric and 
higher pressures in hydrogen and nitrogen free of electron-attaching impurities. This kind of 
tube has advantages where the current demands upon a rectifier are small. 





N the course of some experimental studies of 

corona discharge, methods were found for 
producing and maintaining conditions for free 
electron discharge in gases at pressures of the 
order of magnitude of atmospheric and higher. 
The very large difference in mobility between free 
electrons and positive ions is illustrated in Fig. 1, 
in which are shown the negative and positive 
corona current-voltage characteristics for a 1-mil 
tungsten wire-end discharging in hydrogen at 
atmospheric pressure towards a plane electrode 
at a distance of 3 cm. Characteristics for wire 
sizes between 0.4 and 3 mil, and for distances to 
the plane electrode between 0.3 cm and 10 cm 
were measured and are summarized in an earlier 
report.! 
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Fic. 1. Current-voltage characteristics for negative 
and positive corona. 2 


Now at the National Bureau of Standards. This work 
was performed in 1939-1941 at the Electronic Research 
Corporation, Newark, Ohio. 

' Willard H. Bennett, Phys. Rev. 58, 992-997 (1940), 
see Figs. 2 and 3, p. 994. The abscissae in those figures 


* 


VOLUME 18, MAY, 1947 


As will be seen from the figure, the very large 
difference between the mobilities of electrons and 
positive ions is quite unusual for corona dis- 
charges. These can only be attained by using 
means for eliminating electron-attaching com- 
ponents from the gas in the tube to a degree 
estimated as less than one part per hundred 
million and keeping the gas free of these com- 
ponents to this degree. 

Cold cathode rectifiers were developed taking 
advantage of this large difference in mobilities. It 
is believed that these rectifiers have peculiar 
advantages over conventional types of rectifiers 
in many applications, particularly where the 








a 


Fic. 2. Representative types of tube. (The dimension 64 is 
in inches.) 





were incorrectly labelled and should have been ‘‘kilovolts 
difference” instead of “‘log voltage difference.” 
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Fic. 3. A representative type of rectifier using ten tubes 
of the type shown in Fig. 2b. The condensers are 0.5 
microfarad each and rated at 20 kilovolts d.c. each. The 
multiplier can be a spiral ink line resistor on a Bakelite or 
formica cylinder with 100 megohms per stage. 


current requirements are small, and that a 
description of this kind of rectifiers would be of 
interest. 


GENERAL DESCRIPTION 


Some representative types of tube are illus- 
trated in Fig. 2. In these tubes, the discharge 
electrodes are one- or two-mil tungsten wire- 
ends. The collecting electrode in each is an 
Aquadag coating on the glass wall shown as the 
shaded area in the broken sections. 

The tubes in Figs. 2a and 2b are filled with 
hydrogen at atmospheric pressure. The tube in 
Fig. 2c is filed at eight atmospheres with 30 
percent hydrogen and 70 percent nitrogen. 

Because there is no filament heating require- 
ment with a tube of this kind, rectifiers built 
using these tubes are simpler, lighter, and more 
compact, than rectifiers using conventional types 
of tubes. As an illustration, Fig. 3 shows a 
rectifier supplying 100 kilovolts using a 10 
kilovolt transformer in a Cockcroft-Walton cir- 
cuit. In this kind of circuit, if conventional types 
of tubes were used, the need for filament heating 
current at each stage would make necessary a 
filament transformer with high voltage insulation 
between primary and secondary for each stage. 
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Such a requirement would make necessary a 
bulkier and more complicated equipment than 
that shown in Fig. 3. 

There is no requirement for filament heat-up 
time and no devices are needed for tube pro- 
tection on this account. 

Rectifiers can readily be designed: with clips 
for holding the tubes so that simply reversing the 
tubes in the clips reverses the output polarity of 
the rectifier. 

Although these tubes do not waste electric 
power in heating a filament, they suffer from loss 
of efficiency from another cause. Under normal 
conditions of operation, there is a small positive- 
ion current on the reverse half-cycle. This effect 
together with the fact that even the high free 
electron mobilities are not sufficient to give as 
high current densities as are obtainable with 
electrons in vacuum, has served to hold these 
tubes to an efficiency of about 20 percent to 50 
percent which is lower than most conventional 
types of rectifier unless no greater output cur- 
rents are required than can be adequately sup- 
plied by these cold cathode tubes. 

In Figs. 4 and 5 are shown some efficiency data 
taken on a rectifier using the type of tube shown 
in Fig. 2c. In Fig. 4, the transformer supplies 10 
kilovolts a.c., and in Fig. 5, the supply is 20 
kilovolts a.c. Observations were obtained for 
hydrogen (shown with full lines) and a mixture of 
30 percent hydrogen and 70 percent nitrogen 
(shown with dashed lines). War work prevented 
completing these data, but from some observa- 
tions in a single wire experimental tube, a good 
estimate was obtained for 20 percent deuterium 
and 80 percent nitrogen (shown with the dotted 
lines). From mobility data* it is to be expected 
that the efficiency would be still higher for a 
mixture of 5 percent deuterium and 95 percent 
nitrogen. For comparison purposes, representa- 
tive curves for a vacuum-tube rectifier are 
included. 


PREPARATION OF A TUBE 
As discharge electrodes, wire-ends. with any 
diameter between about 0.4 mil and 3 mils seem 


2 Additional engineering data is given in Trans. A.I.E.E. 
64, 164-7 (1945); and Electronics Industries, p. 106 
(November, 1945). 

?W. H. Bennett and L. H. Thomas, Phys. Rev. 62, 
41-47 (1942). 
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to be about equally effective. The reason for this 
is that in the initial conditioning of the tube, 
there is an etching of the wire-ends to a sharp 
conical tip whose radius of curvature is of the 
order of magnitude of 0.2 mil.‘ In tubes con- 
taining six or more wire-ends, the maximum 
number that can usefully be grouped with bases 
adjoining is six, and even with six wires in a 
group, the maximum current per wire-end that 
can be drawn without arc-over is less than six 
times the current that could be drawn from ‘a 
single wire, the other electrode spacings being the 
same but the other five wires removed. The dis- 
charge spreads laterally from a wire-tip and 
adjacent discharges interfere with each other to 
some extent.® 

It is important that the electric field diverge 
all the way to the collecting electrode. Aquadag 
or metal collecting electrodes can be used. 


After filling the tube with gas and sealing the 


tube, a.c. voltage is applied directly. The result- 
ing discharge eliminates the electron-attaching 
impurities and sharpens the wire-ends, slowly at 
first. The current rises quite rapidly during the 
last stages of the conditioning, however, and 
unless a circuit breaker is used in the a.c. supply, 
the tube may be greatly over-heated. 

This does not complete the conditioning. If the 
tube is allowed to stand, it will probably decondi- 
tion to some extent after a day, that is, the 
electron mobilities will have decreased by a large 
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Fic. 4. Efficiency with 10 kv a.c., applied. 


‘Willard H. Bennett, Phys. Rev. 61, 55 (1942). See also 


reference 1. 
5 Willard H. Bennett, Phvs. Rev. 61, 53-56 (1942). 
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amount. This is probably caused by the liberation 
into the gas of some electron-attaching impurities 
from the walls of the tube and from the electrodes. 
For this reason, immediately following the first 
run-up, the tube is put on a reduced a.c. voltage 
sufficient to maintain a current of about 100 
microamperes per wire. After about three hours, 
the tube is permanently conditioned. 


OPERATION 


One of the convenient features of this kind of 
tube is that the full rated voltage can be thrown 
on without any preliminaries. 

There is a definite output voltage at which a 
tube arcs over and the maximum d.c. output 
rating of the tube can be given as about 10 
percent under that value. The tube never “‘arcs- 
back”’ at unexpectedly low voltages as some of 
the conventional mercury vapor tubes do 
occasionally. 

There is no observable consumption of dis- 
charge electrodes or other materials in these 
tubes and there is no known reason for these 
tubes to have other than an indefinitely long life 
when operated within the ranges of rated current 
and voltage. 

Surges im the high voltage a.c. supply, caused 
by opening or closing the primary circuit, cannot 
puncture the tube with cold emission because the 
tube is filled with gas. 

The small positive-ion current on the inverse 
half-cycle serves a useful purpose. When the 
output of the rectifier is being smoothed with a 
condenser, resetting the input a.c. voltage to a 
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Fic. 5. Efficiency with 20 kv a.c., applied. 
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w . 
lower value results in the condenser voltage’ 
dropping correspondingly because of the positive- 
ion current in the tube. 


As long as the rated inverse voltage is not 
exceeded with these tubes, no damage results 
from drawing more than the rated current through 
the tubes momentarily. The only damage that 
eventually occurs is caused by over-heating the 
tube as a whole and this takes at least several 
minutes for all types of tubes. The largest tubes, 
containing part nitrogen, decondition at currents 
of about 0.2 milliampere per wire because of a 
slow formation of ammonia. Operation at 0.06 
milliampere per wire slowly disassociates the 
ammonia and reconditions the tube. 

At the higher values of current, ionization by 
collision occurs between the sheath and the col- 
lector in a tube and a high pressure glow dis- 
charge makes its appearance. This glow discharge 
has a positive resistance characteristic for cur- 
rents up to a higher order of magnitude than the 
limiting current for the corona discharge, and is 
not to be confused with normal high pressure 
arcs. This is not a disruptive kind of arc (that is, 
it is not abruptly initiated with a streamer and a 
spark) because of the absence of electron- 
attaching gaseous components. ® 

In tubes containing part nitrogen, at currents 
of the order of 4 milliampere per wire-end, the 
wire-ends become heated, and at currents of be- 
tween 1 and 5 milliamperes per wire-end, the 
~ Willard H. Bennett, J. App. Phys. 13, 199-200 (1942). 


William H. Otto and Willard H. Bennett, J. Chem. Phys. 
8, 899-903 (1940). 


wire-ends attain temperatures as high as 2000°C. 
The energy necessary to produce these tempera- 
tures cannot be accounted for by the impact 
energies of positive ions striking the wire-ends 
following their last free paths. Such heating does 
not occur in discharges in pure hydrogen at the 
same currents.’ It appears that either an active 
state of nitrogen, or else an excessive concentra- 
tion of atomic hydrogen caused by the presence 
of nitrogen, is produced in the ionization sheath 
and is responsible for this excessive heating 
which occurs only in part nitrogen. 

When the wire-ends heat to more than about 
1000°C, a reduction in sheath drop in potential! 
occurs both for positive and negative discharge 
because of the thermal reduction in gas density 
immediately adjacent to the wire-end. There is a 
corresponding reduction in arc-over voltage for 
the tube on the inverse half-cycle. This effect is 
reduced by using the finest wires (0.4 mil) but is 
not entirely overcome by this means. 


CONCLUSION 


It is not believed that the full potentialities of 
this kind of discharge in a rectifier have been 
realized. Further studies of the diffuse kinds of 
arcs which occur in these gases should prove 
profitable. It is thought that a wider appreciation 
of the advantages of a cold cathode gas-filled 
rectifier may bring an increased interest in 
studies of the properties and methods for control 
of these arcs. 


~ 7 Willard H. Bennett, Phys. Rev. 62, 369-371 (1942). 





The Hydrodynamic Lubrication of Finite Sliders 


CHARLES P. BOEGLI 
Springfield, Ohio 
(Received December 17, 1946) 


Two approximations are made in the solution of Reynolds’ lubrication equation for the case 
of a finite slider. These approximations iead to a series of equations that are easy and rapid to 
use for flat sliders, and which are also applicable to curved sliders. Examples of their use are 
presented, and a number of calculations are made to determine the range of slider proportions 

; to which they may be applied with sufficient accuracy. The method is of special value for 
wide sliders which are difficult or impossible to analyze by other methods. 


INTRODUCTION 


HE solution of the basic Reynolds differ- 
ential lubrication equation for the case of 
a slider bearing of finite width has always been 
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beset with great difficulties. Reynolds! did not 
succeed in finding a general solution to the 


10. Reynolds, “On the theory of lubrication and its 
application to Mr. Beauchamp Tower's experiments,”’ 
Phil. Trans. Roy. Soc. 177, part 1, 157-234 (1886). 
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equation, and it was not until nineteen years 
after his presentation that Michell? provided a 
satisfactory attack to the problem. 

Michell’s method involves no assumptions 
that affect the accuracy of the solution for the 
specified conditions. It is, however, applicable 
only to sliders with flat surfaces, and even for 
such simple cases the solution is so tedious that 
only a very small amount of data has yet been 
computed. R. O. Boswall, working with an 
assumption credited to Stodola, succeeded in 
developing an approximate method* that gives 
results negligibly different from those of the exact 
method in cases where comparisons are possible. 
Here again the analytical solution can only be 
made for flat sliders; for other shapes resort is 
had to a graphical method that actually consists 
in estimating by eye the shape of one of the 
pressure distribution curves. 

More recently, Muskat, Morgan, and Meres 
published another solution to the original differ- 
ential equation. They presented their solution 
as a series of graphs that serve better for the 
solution of problems of operation than problems 
of design. Again, only bearings with flat surfaces 
are considered. 

There is great need for a set of equations that 
will provide, with a minimum of effort, a good 
approximation to the actual behavior of a wider 
variety of sliders than can now be handled. The 
need for a very high degree of accuracy is, for 
several reasons, not great. One reason is that it is 
never possible to predict closely the behavior of 
the viscosity of the lubricant as it passes through 
the bearing ; the assumption of constant viscosity 
is generally made, and the practical value of the 
solution is limited by this approximation. Also, 
it would be very difficult actually to construct 
a bearing in exact conformity with a specified 
design ; such a small difference, for example, as a 
slight rounding of the leading edge, would have 
considerable effect upon the performance. There- 
fore, a solution that predicts to a fair degree of 
accuracy the performance of a great many types 

2 A. G. M. Michell, ‘‘The lubrication of plane surfaces,” 
Zeits. f. math. Physik 50, 1073-1099 (1904). 

> R. O. Boswall, The Theory of Film Lubrication (Long- 
1 i Green, and Company, New York, 1928), pp. 141- 


* Muskat, Morgan, and Meres, ‘“‘Studies in lubrication. 
VII. The lubrication of plane sliders of finite width,” J. 
App. Phys. 11, 208-219 (March 1940). 
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of slider bearings would be at least as valuable 
as one that predicts very closely the performance 
of only a narrow range. 

Such a solution is provided through two simple 
assumptions. The first is that used by Boswall 
in his derivation of an approximate method; 
namely, that the pressure functions along the 
length and width of the slider are independent. 
The second is that the pressure function along 
the length of the slider (that is, the dimension 
in the direction of motion) is the same as that 
for an infinite slider. 

This paper presents analyses based upon these 
two assumptions; an analytical solution for the 
case of the flat slider and an arithmetical method 
for curved sliders. Since with any approximate 
method some estimate of the accuracy and range 
of applicability is desirable, the report also 
presents a comparison of the results obtainable 
with the new method to any available data 
found in other manners. 


ANALYSIS 


The general differential equation of Reynolds 
serves as the starting point for this analysis of 
slider lubrication. 


d/10P d/10P d 
—(- —i*) +—(- —k*) =6—Uh. (1) 
Ox\u Ox dz\p dz dx 


In this equation yu, h, U, and P are the oil vis- 
cosity, the oil-film thickness, the slider velocity, 
and the oil pressure at any point on the slider 
surface. If the substitutions are made: 


x=xXot+Lnx,, 
2= Bai, 
h —_ ho, 
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where L and B are the length and width of the 
slider, 4o represents the oil-film thickness at the 
outlet edge, and £ is a function of x; (see Fig. 1) ; 
and if viscosity is assumed constant, the equation 


becomes 


a aP a°P dg 
G [eave =A (2) 
1 


Ox Ox; 


In Eq. (2), N is the ratio of slider length to slider 
width (L/B). and A is the product 64¢UL/h¢ 
containing all the dimensional terms. The new 
variables x; and z; are the fractional distances 
along the length and width of the slider, and 
range in value from 0 to 1.0. 

It will now. be assumed, following the method 
in reference 3, page 141, that the pressure at any 
point can be expressed as the product of two 
functions 


P=A-f(x1)-f(s1). (3) 


This product is substituted for P in Eq. (2), 
yielding 


a] df ( x1) d*f(z1) dp 
. co — |+ N7B*f(x1)— - (4) 


Ox ax, dz," dx, 


To find the form of f(z:) Eq. (4) is solved along 
the line where f(x,) is a maximum; that is, where 
df(x,)/dx, is zero. This yields 


1 a] 1 =| 
d*f(z1) f(x.) dx," Im B*f(x1) dx, Im 
cae of f 











dz? N? - Nn 

(5) 
in which x, indicates the value of x; for which 
f(xi) is a maximum. This equation can be 
simplified by considering that since f(x,) is the 
same as that for an infinite slider, this being the 
. second basic assumption, this function was de- 
rived from the differential equation 


d dP dp 
£6). 
dx, dx, dx, 
which, upon substitution of P=A -f(x,) and per- 
formance of the indicated differentiation, gives 
df (x1) af (x1) dp 
36%*°——_+- 6° —__ = - 


dx, dx," dx, 
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Where f(x1) is a maximum, this equation shows 


that 
d2f(x;) 1 dg 
dx, Im B* dx, Im 


and consequently Eq. (5) is reduced to 
d*f(z:) MM M 

—+-—_f(s;) = — 

dz;> —- N? N? 


where J/ is equal to the expression on either 
side of Eq. (6) divided by f(x»). The solution is 


e*—1 
fi) =1-( ~)ee mie —a(l—s)), = (7) 
e 


rn—@ n 


where 


Up to this point the treatment has been 
general, and applicable to both flat and curved 
sliders. In the case of flat sliders, 





B= 1 + bx, 
and 
(b+2)! 
M = —}——— 
(b+1)? 
while f(x1) is simply5 
bx\(x1—1) 
f(x1) =— — (8) 
(b+2)(bx,+1)? 


For curved sliders, an arithmetic or graphical 
evaluation of f(x:) is usually preferred to an 
analytical solution. Such methods result in a 
series of tabulated values of the function for 
various values of x;, and it is necessary to find 
Xm and f(x») from this table. An excellent method 
of finding the required values is presented in 
Empirical Formulas.® To evaluate n, the value of 
V must be known; it can be calculated from the 
expression on the right side of Eq. (6), which is 
to be preferred to the quantity on the left, 
because the former does not require the finding 
of second derivatives from graphical data. 

The important equations and constants are 


5 Reference 3, page 44. 
6 T. R. Running, Empirical Formulas (John Wiley and 
Sons, Inc., New York, 1917), pp. 100-113. 
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Fic. 2. 
summarized below: 
P=A -f(x1)-f(21) (3) 
e"—] 

f(a) = -( Je msif-e-i-s)) = (7) 

e"—e-” 

bxy(x1—1) 





f(x) =— ties, a for flat sliders only 
(6+2)(bx,+1)? 
(8) 
in which 
6uUL 
A =—— 
ho? 
M 
n? = ——— 
N? 
- 1 d*f(x1) ~ 1 dp 
meee 
F(x) dx) Im B*f(x1) dx, Im 
(b+2)! 
= — }———_ for flat sliders only 
(b+1)? 
N=L/B. 


EXAMPLES OF FLAT-SLIDER ANALYSIS 
To illustrate briefly the use of the equations, 
let it be desired to analyze the lubrication of a 
slider for which N=1.0 and 6=1.5. For this 
slider, 17J=—12.0 and mn is thus found to be 
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3.46. The pressure functions are then 
1.5x1(x1— 1) 
f (x1) = 
(3.5) (1.54141)? 
f (21) =] —().97(e 3.46814. ¢-3.46(1—s1)) | 





Graphs of these two functions can be drawn, and 
by finding curves of constant f(x,)-f(z:) the 
isobars over the slider surface can be located. 
This has been done in Fig. 2. The numbers on 
the isobars are values of f(x1)-f(z:). The actual 
pressure at each isobar is, of course, this value of 
f(x1)-f(z1) multiplied by the quantity A. 

For the case of wide sliders, where N becomes 
small, the analysis by means of the rigorous 
equations, or even by the approximate treat- 
ment in reference 3 becomes more difficult, but 
that with the new equations becomes simpler. 
For a slider where N=0.01 and b=1.5, M is 
still —12.0 but » is now 346. Then, 


f (21) _ 1 a (e~ 34621 4. e—346(1—21)) | 


It is observed that the value of this function, 
although zero at 2,=0 and 1.0, rises very 
rapidly to 1 between these endpoints; at 2; =0.02 
or 0.98, for example, the value of f(z:) is 0.999. 
The function f(z;) may therefore be considered 
as a correction to be applied to the ends of a 
slider which would otherwise be considered in- 
finite. Under such conditions, the value of f(2:) 
should be directly dependent upon 2;/N, the 
actual width represented by 2; that this is 
actually the case is illustrated by finding the 
point at which f(z;) =0.999 for a slider for which 
N=0.0001 and b=1.5. This value turns out to 
be z,=0.0002, which is 0.01 times that for a 
slider of N=0.01. Since, however, this second 
slider is 100 times as wide, the actual distance 
from the ends of both sliders to the points at 
which f(z;)=0.999 is the same. This type of 
relationship in fact exists up to the point where 
the quantity 


(e” — 1 ) (e” —e-") 
becomes significantly different from 1.0. 
EXAMPLE OF CURVED-SLIDER ANALYSIS 


For purposes of illustration, a slider will be 
analyzed for which B= (2.50)** and N=0.5. The 
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customary approximate method leads to a series 
of values for f(x1): 


x1 B f (x1) 

0.0 1.000 0.0000 
0.1 1.096 — 0.0273 
0.2 1.201 — 0.0414 
0.3 1.316 — 0.0466 
0.4 1.443 — 0.0459 
0.5 1.582 — 0.0414 
0.6 1.733 — 0.0349 
0.7 1.899 —0.0273 
0.8 2.081 — 0.0192 
0.9 2.281 —0.0112 
1.0 2.500 0.0000 


By the use of Lagrange’s interpolation formula, 


TABLE I. 
b=1.0 b=2.0 
W/W. W/W. Error W/Wa W/W. Error 

N (new) (actual) % (new) (actual) Qq 
0 1.000 1.000 0 1.000 1.000 0 
0.5 0.687 0.695 —1.2 0.735 0.714 2.9 
1.0 0.421 0.434 —3.0 0.493 0.460 7.2 
1.5 0.259 0.279 —7.2 0.322 0.297 8.4 
2.0 0.168 0.185 —9,2 0.219 0.208 5.3 
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Xm is found to be 0.338 and f(x,,) = —0.0470. At 
this point B=1.363 and dB/dx,;=1.249 


1.249 
M =—————————__- = — 10.50. 
(—0.0470) (1.363)? 


The value of n is then 
n = (10.50)!/0.5 = 6.482. 


All the necessary constants have now been calcu- 
lated. The shapes of the isobars over the surface 
are illustrated in Fig. 3. The value marked on 
each isobar gives, when’ multiplied by A, the 
magnitude of the pressure at that isobar. 


APPLICABILITY OF EQUATIONS 


Although the new equations are especially 
useful in analyzing cases that are difficult or 
impossible to solve by the more exact and in- 
volved methods, the only available accurate 
data are for flat sliders, and these must conse- 
quently serve as a basis for estimating the 
accuracy of the solutions obtained by the new 
method. The comparison thus made is presumed 
to hold also for curved sliders. 


W/W. Ratio 


By W/W, ratio is designated the ratio of the 
total load carried by the slider to that which 
would be carried by an equal length of an 
infinitely wide slider. Since f(x,) in the new 
equations is the same as that for an infinite 
slider, the W/W, ratio is merely the integral 
of f(z:) (Eq. (7)), which is 


W 2 (1-—e—")? 
— =|-—~- ——_. (12) 
W. n (1—e-?") 

Values computed from this equation are com- 
pared to those computed by Michell’s method 
for b=1.0 and 2.0 in reference 3 
(Table 1). 


, page 137 


TABLE II. 
b=2.0 N D (new) D (actual) Error % 
0 — 0.0246 — 0.0246 0 
0.5 —0.0181 —0.0175 3.4 
1.0 —0.0122 —0.0113 8.0 
‘5 — 0.0079 — 0.0073 8.2 
2.0 —0.0054 —0.0051 5.9 
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It will be shown later that the error in the 
W/W. ratio is of the same magnitude as that in 
many other derived characteristics of the slider 
bearing. Figure 4 compares values of W/W, 
obtained with the new equations to those found 
by several investigators with various methods, 
for N=0.5, 1.0, and 2.0. The values due to 
Kingsbury are from “Optimum Conditions for 
Journal Bearings’”’ and those credited to Muskat, 
Morgan, and Meres were actually found by a 
series of graphical calculations with the charts 
in Lubrication,’ pages 162-66. The latter cannot 
therefore be considered extremely accurate, but 
they extend over a greater range than the other 
published values. 

It can be seen that the ratio given by the new 
equations is in good agreement with those com- 
puted by other methods. 


Total Load 


The total load carried by the slider is given by 
the expression 


W=ALBD, 


following the notation of Boswall. The new 
method yields, for the value of D, 


p= f ssdds fi sod, 


because the variables x; and 2; are independent. 
Comparison of the values of D from the new and 
from the Michell solution is sufficient because 
the factors A, L, and B are the same for both 
methods. The accurate data are again obtained 
from reference 3, page 137 (Table II). It could 
have been expected that the error in total load 
be somewhat greater than that in W/W, ratio, 











TABLE III. 
b=2.0 N xi (new) x1 (actual) Error % 
0 0.392 0.392 0 
0.5 0.392 0.386 1.5 
1.0 0.392 0.374 4.9 
1.5 0.392 0.359 10.5 
2.0 0.392 0.341 12.6 





7 Albert Kingsbury, ‘Optimum conditions for journal 
bearings,”” Trans. ASME RP-54-7, 54, 123-148 (1932). 

8A. E. Norton, Lubrication (McGraw-Hill Book Com- 
pany, Inc. New York and London, 1942), pp. 162-166. 
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because here the variation of f(x,) from the 
actual form is also taken into account. Never- 
theless, the accuracy is still surprisingly good. 


Center of Pressure (Z;) 


The center of pressure is the value of x; at 
which the integrated effect of the pressure is 
directed. It is thus the line along which the 
slider would be pivoted. This depends, in the 
new method, only upon f(x.) which keeps the 
same form regardless of N. The value of Z, is 
therefore a constant depending only upon 6, and 
is equal to #, for an infinite slider. It is, of course, 
known that the value of Z; is actually dependent 
upon N, but the error introduced by assuming it 
constant remains to be found. The actual values 
are obtained from the same source previously 
quoted (Table III). The error in &, is still not 
excessively great, and the assumption of the 
constant value would lead to no serious difficulty 
in most cases. 




















TABLE IV. 
b=2.0 N H (new) H (actual) Error % 
0 1.82 1.82 0 
0.5 2.00 2.02 —1.0 
1.0 2.30 2.38 —3.4 
1.5 2.74 2.83 —3.2 
2.0 3.23 $.36 —2.7 
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TABLE V. 


b=2.0 N Vo (new) Vo (actual) Error % 
0 0.288 0.288 0 
0.5 0.225 0.227 —0.9 
1.0 0.169 0.172 —1.8 
1.5 0.126 0.130 —3.1 
2.0 0.100 0.115 — 14.0 

b=2.0 N Vi (new Vi (actual Error % 
0 0.288 0.288 0 
0.5 0.313 0.313 0 
1.0 0.306 0.299 ye 
1.5 0.274 0.262 4.7 
2.0 0.240 0.245 2.0 


Frictional Force 


In reference 3, page 145, the frictional resist- 
ance is given by the equation 


woU\? 
R=W H, 

W, 
where /H is an expression involving integrals of 
f(x.) and f(z,). As in the case of total load, a 
comparison of values of 7 is sufficient, since the 
remaining terms are identical for both methods 
(Table IV). Here, the error is very small indeed. 


Quantity of Lubricant 


The following expressions are given for the 
quantity of lubricant in reference 3: 


; B,° : 1 ; 
| i=4(6D)! Bi-- G: f f(zi)dz, 
€ 0 


; ea 
Gi= meomes ’ 
dx, r=1 
1 
Vo= 46D) 1-Go f fiends, | 


- —"] 
Go= ° 
dx, z1=0 


V, and Vo are the inlet and outlet quantities, 
respectively. Substitution of the proper values 
into these equations yields the results given in 


where 


where 
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TABLE VI. 
b=2.0 N W/W. D x1 H Vi Vo 
0 0 0 0 0 0 0 
0.5 2.9 3.4 1.5 —1.0 0 —0.9 
1.0 7.2 8.0 49 —3.4 a3 —1.8 
1.5 8.4 8.2 10.5 —3.2 4.7 —3.1 
2.0 5.3 5.9 12.6 —2.7 2.0 — 14.0 


Table V. The error in Vo is too large at N=2.0 
but is certainly sufficiently small at N=1.5. The 
error in V,; is small enough over the entire 
considered range of N. 


Estimation of Error 


It was noted previously that the value of 
W/W. is very easily calculated by the new 
method, and it would be desirable to use the 
error in this quantity as a criterion of the 
accuracy of the method. For purposes of com- 
parison, then, the errors are tabulated in Table 
VI. It is considered that the W/W, ratio forms 
a fair, if not slightly pessimistic, criterion of the 
accuracy of D, H, and V,, while for ; and Vo the 
error is likely to be greater in magnitude than 
that estimated from W/W, when N becomes 
large. For the great majority of sliders, for which 
N&1.0, the W/W, criterion is satisfactory. 

Since Fig. 4 shows good agreement between 
actual and computed W/W, ratios even at 
N=2.0, it can be concluded that the new equa- 
tions hold with good accuracy for all sliders for 
which N is 1.0 or less, and they provide an 
excellent basis for estimating the behavior of 
even narrower sliders. 


ACKNOWLEDGMENTS 

The original work on this project was carried 
out at the National Advisory Committee for 
Aeronautics, Aircraft Engine Research Labora- 
tory, Cleveland, Ohio, and was published as 
ARR E4122. The present report constitutes a 
refinement in treatment and a considerable ex- 
pansion of that work. 

The assistance of Miss Thelma Bresser in 
reading the proofs is gratefully acknowledged. 


JOURNAL OF APPLIED PHYSICS 











ie 


Letters to the Editor 











Measurements of Base Pressure on a Missile in 
Free Flight* 


F. K. HILL AND R, A. ALPHER 


Applied Physics Laboratory, The Johns Hopkins University, 
Silver Spring, Maryland 


February 24, 1947 


N a series of experiments now under way, missiles of 

simple geometry (see Fig. 1), differing only in length, 
are rocket-launched to supersonic velocities; telemetering 
installed in the missiles transmits to a ground receiver the 
aerodynamic data measured during flight. Among the 
significant data being obtained are base pressures; such 
data are important in the extrapolation of wind-tunnel 
studies on models to full-scale performance of missiles, 
where large Reynolds numbers are encountered under 
ordinary atmospheric conditions. 

The missiles have a 14-caliber ogival nose joining tan- 
gentially to a cylindrical body of 4-inch diameter, the base 
of the cylinder being flat and closed. Four fins of bi-convex 
section are added to give flight stability. The yaw of the 
missile (in the plane:including the missile axis and base 
pressure measurement port) with respect to the flight path 
is determined by differential pressure measurements near 
the nose tip. These differential pressure measurements 
may be correlated with yaw by means of laboratory experi- 
ments in an open supersonic jet or in wind tunnels, as well 
as by calculations based on linearized supersonic flow 
theory. It is found in these experiments that the yaw 
angles encountered are of the order of 1° or less and that 
no measurable effects caused by angle of attack are observed 
on base pressure data here reported. Missile velocities are 
obtained both by an optical tracking system and by 
stagnation pressure measurement on the missile with a 
Pitot tube. 

Base pressures were measured at a 0.040” port located 4” 
from the center line of the missile. Observed base pressures 
were converted to pressure coefficients by dividing them 
by the dynamic head (namely, 4Xdensity Xvelocity?); 
corrections caused by change in altitude during flight were 
made. Coefficients thus computed are shown (Fig. 2) for 
missiles of the first two lengths of the series to be studied. 
It is to be noted that the coefficients show at most a small 
variation with Mach number in the supersonic range 
covered by the tests, where the value is about 0.3. This 
result is significantly different from previously used theo- 
retical values.'! While present information indicates that 
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Fic. 1. Scale effect test missile. Telemetering to record: (1) Pitot tube 
stagnation pressure, (2) differential pressure at nose ports, (3) base 


pressure at port $” from center. All pressure ports in the same plane 
through missile axis. 
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Mach Number 


Fic. 2. The probable errors indicated are based upon fairly well-known 
inaccuracies in the determination of the dynamic head and static pres- 
sure throughout flight, and an assumed inaccuracy in the telemetering 
records based on laboratory tests. All errors are believed to be random. 


pressure is essentially independent of position on the base 
of such missiles, an experimental check of this will be made 
in free flight. 

Similar missiles with lengths up to 120”, as well as a 
series differing only,in diameter, and a series of homologous 
models, are to be studied, with maximum velocities corre- 
sponding to Mach number 3. 

This study is part of a general investigation of scale 
effects which is to be reported in detail after completion of 
all the tests. It provides the first direct measure of base 
pressure on such high speed missiles. 

* Work being performed by a group under contract NOrd 7386 with 


the Bureau of Ordnance, U. S. Navy 
10.2 given by von Karman and Moore, Trans. ASME 54, 303 (1932). 





Direction of Maximum Crystal Elongation 
during Metal Cutting 


G. H. TowNEND 
Metropolitan-Vickers Electrical Company, Ltd., Manchester, England 
March 17, 1947 


N an interesting paper by Dr. M. E. Merchant on 

““Mechanics of the Metal Cutting Process,”"! an error 
was made in the expression connecting shear angle, rake 
angle, and direction of maximum crystal elongation (Eq. 
(1), page 269). I have written to Dr. Merchant and he 
confirms the new expression derived below. 


Let a=rake angle, 
¢ = the angle between the shear plane and the surface 
of the workpiece. 
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Consider a section of the workpiece ABCD such that 
Z CBD =90°. Then after shearing this section will become 
the parallelogram BFEC. The dimension perpendicular to 
the shear plane must be unchanged during shearing. 


BD=BH=a. 
Consider a line CG in the chip making an angle ¥ with the 
plane of shear. Then the corresponding line in the work- 
piece is DJ. The condition that CG shall be the direction 
of maximum crystal elongation is that CG/DJ is a 
maximum. 
CG =a cosecy, 
DJ = \a*+(a cot¢—a coty +a tan(¢—a)*}}. 
The condition is that 
{sin’*y + (cote siny —cosy +tan(@—a) siny )?}-4 
is a maximum. That is, 
sin*y + {siny (cote +tan(¢@ — a) — cosy }? 
is a minimum. Let 
cote + tan(@—a) =e, 
é sin2y —2¢ cos2y =0. 
If «+0 
2 cot2y =e, 
2 cot2y =cot¢+tan(o—a). 


Alternatively it is easily shown that a circle on the 
workpiece becomes an ellipse on the chip and that the 
principal axes of the ellipse make an angle y with the shear 
plane, where 

2 cot2y =cot¢+tan(o—a). 


1 J. App. Phys. 16, 267 (1945). 
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Printing Electronic Circuits 


According to a recent release from the National Bureau 
of Standards prepared by Cledo Brunetti and Philip J. 
Franklin, new techniques, developed by the National 
Bureau of Standards in cooperation with electronic manu- 
facturers, for electronic devices promise to simplify pro- 
duction methods and to accelerate the manufacture of 
cheaper, more compact, and more uniform electronic equip- 
ment. These techniques, representing the perfection of a 
method of printing wiring and circuit components on an 
insulated surface, is one of the new processes and tech- 

‘niques evolved from research and development during 
the war. 

The methods are applicable in the design of numerous 
industrial and commercial devices where extreme rugged- 
ness and small size are imperative. By means of the printing 
process the electronic circuit of a conventional radio— 
now characterized by a maze of wires, resistors, inductors, 
and condensers—may be replaced by a compact circuit 
printed on a small flat plate, eliminating much of the third 
dimension and making it simple to check and repair. 
Other specific applications lie in the fields of pocket radios, 
personal radio telephones, miniature hearing aids, meteoro- 
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logical instruments, and various miniature electronic con- 
trol circuits. 

Present investigations at the National Bureau of Stand- 
ards of the various means available for reducing the size 
of electronic assemblies and for mechanization of chassis 
wiring reveal a number of feasible and practical electronic 
printing methods. Variations of the stenciled screen-print- 
ing process include hand printing or spraying of the com- 
ponents and leads on suitable surfaces. Another method 
involves the use of photography and an abrasive blast. A 
third is an adaptation of the chemical methods of silvering 
mirrors. Mechanical processes include metal spraying, 
metal sputtering or evaporation, the use of foils, and stamp- 
ing. As an example, loop antennas are now stamped out 
in a single operation by a variation of the printed wiring 
techniques and have shown improved performance over 
wire-wound loops in addition to reduced cost of production. 
The principal effect of these methods is to reduce electronic 
circuit wiring to two dimensions. The effect is enhanced 
where it is possible to employ sub-miniature tubes and 
compact associated components. 


New Appointments 


Leonard I. Schiff, theoretical physicist who worked on 
the Los Alamos atom bomb project, has joined the faculty 
of Stanford University as associate professor of physics. 

H. B. G. Casimir, co-director of the Philips Research 
Laboratories, Eindhoven, Holland, is delivering a series 
of invitation lectures at Johns Hopkins University on the 
properties of matter at low temperatures and on problems 
in quantum electrodynamics. 

Leo L. Beranek, formerly director of the electro-acoustic 
and systems research laboratories of Harvard University, 
has been appointed associate professor of communications 
engineering in the Department of Electrical Engineering 
of the Massachusetts Institute of Technology. 

Alfred J. Reis, Austrian-born scientist, was recently 
appointed a research specialist in the Engineering Experi- 
ment Station at Rutgers University. He will conduct a 
series of studies in the engineering properties of materials 
and in powder metallurgy. 

John A. Hipple has been named chief of the Atomic 
Physics Section at the National Bureau of Standards. 

John L. Abbott recently joined North American Philips 
Company, Inc., as application engineer in the industrial 
x-ray division. 


Awards 


James I. Hoffman of the National Bureau of Standards 
has been awarded the 1946 Hillebrand Prize for significant 
contributions to chemical science. 

The John Price Wetherill and the Howard N. Potts 
medals of The Franklin Institute of Philadelphia were 
awarded at ceremonies on April 16. Kenneth S. M. David- 
son of Stevens Institute of Technology received the 
Wetherill medal for his contributions to the improved 
design of ships through experiments on ship hull models. 
Vladimir K. Zworykin, director of electronic research of 
the Radio Corporation of America, and Robert H. Kent, 
associate director of the ballistics research laboratories, 
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Aberdeen Proving Ground, were each awarded a Potts 
medal. Zworykin received the medal for his inventions of 
the iconoscope and kinescope, which are essential to 
modern commercial television. Kent was honored for con- 
tributions to the science of ballistics and particularly for 
his effective operation of the Aberdeen Wind Tunnel, 
which is doing a large part of the supersonic testing for 
the armed forces. 

The seven winners of Frank B. Jewett fellowships for 
research in the physical sciences for 1947-1948, announced 
by the American Telephone and Telegraph Company, are 
as follows: M. G. Ettlinger of California Institute of 
Technology and Austin, Texas; Wallace D. Hayes of 
California Institute of Technology and Palo Alto, Cali- 
fornia; Paul Olum of Harvard University and Winchester, 
Massachusetts; Aadne Ore of Yale University and Oslo, 
Norway; Alfred Schild of the Carnegie Institute of Tech- 
nology and Pittsburgh, Pennsylvania; Robert L. Scott of 
the University of California and Santa Ana, California; 
and Edwin H. Spanier of the University of Michigan and 
St. Paul, Minnesota. Two of the winners are chemists, two 
physicists, two mathematicians, and the seventh is an aero- 
nautical engineer. 


General Electric-Yale University Joint Research 


At the Dunham Laboratory of Electrical Engineering 
at Yale University, research is being conducted on controls 
for a flight simulator or dummy rocket for duplicating 
flight characteristics of a rocket missile, as part of General 
Electric Company’s long-range rocket research program. 
The project at Yale is under the direction of John L. 
Bower, Assistant Professor of Electrical Engineering. 


Biological Photographic Exhibition 


An international exhibition of biological photography of 
interest to all photographers in the scientific field will be 
held September 10-27, 1947, under the sponsorship of the 
Biological Photographic Association at the Rochester 
Public Library, Rochester, New York. Three judges who 
are prominent in their fields will select all entries. Member- 
ship in the Association is not a requisite for participation 
in the exhibition. Entry blanks for contributors may be 
obtained from Mr. H. Lou Gibson, 343 State Street, 
Rochester 4, New York. The entry fee is $1, and the last 
day for receiving entries is August 1, 1947. 

Information concerning the convention program of the 
Biological Photographic Association to be held September 
10-13 may be obtained from Mr. John V. Butterfield, 
635 St. Paul Street, Rochester 4, New York. 


General Electric Educational Fund Awards 


Research fellowships totaling $20,500 have been awarded 
to fourteen graduate students by the General Electric 
Educational Fund. Eight of the winners received Charles 
A. Coffin Fellowships awarded for research in the fields of 
electricity, physics, and physical chemistry. The other six 
young men obtained Gerard Swope Fellowships for ad- 
vanced study in industrial management, engineering, the 
physical sciences, and other scientific and industrial fields. 
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M.I.T. Summer Session in Applied Mathematics 


Massachusetts Institute of Technology has announced a 
special summer session in applied mathematics from 
August 4 to September 19, 1947. The following graduate 
courses will be offered: 

1. Theory of Plates and Shells. 

AND E. REISSNER. 


2. Theoretical Hydromechanics. By Proressors E. REISSNER AND 
Cc. C. Lam. 

3. Advanced Topics in Applied Mathematics. By Proressors J. L. 
SYNGE, H. REISSNER, W. PRAGER, F. B. HILDEBRAND AND P. D. 
Crovut. 

4. Tensors in Mechanics. By Proressor D. J. StruIk. 


By Proressors F. B. HILDEBRAND 


Further information concerning the program may be 
obtained by writing to Department of Mathematics, 
Massachusetts Institute of Technology, Cambridge 39, 
Massachusetts. 


Centenary of Association of Engineers 


The Association of Engineers, with headquarters at the 
University of Liége in Belgium, celebrates this year its 
100th anniversary. For the occasion it is sponsoring a 
congress and an exhibition to which engineers and industry 
representatives are cordially invited. The congress will 
take place in Liége from August 30 to September 13, and 
is entitled “‘Past, Present and Future of Our Industry.” 
For further information address Association des Ingenieurs, 
12, Quai Paul Van Hoegaerden, Liége, Belgium. 


National Electronics Conference 


The National Electronics Conference, Inc., whose pur- 
pose is to serve as a national forum on electronic develop- 
ments and their application, is sponsored jointly by Illinois 
Institute of Technology, Northwestern University, Ameri- 
can Institute of Electrical Engineers, Institute of Radio 
Engineers, and the University of Illinois, with the Chicago 
Technical Societies Council a cooperating organization. 
Plans are now being formulated for the 1947 National 
Electronics Conference, which will be held on November 3, 
4, and 5 at the Edgewater Beach Hotel, Chicago. 


Course in Statistical Quality Control 


An advanced eight-day course in quality control by 
statistical methods is to be presented June 27 through 
July 5 at Purdue University, Lafayette, Indiana. It is 
intended for those who already have some acquaintance 
with statistical quality control from such a course as the 
War Production Board organized during the war. Anyone 
desiring further information should write to Professor 
I. W. Burr at Purdue. 


Marconi Anniversary 


An International Congress will celebrate the fiftieth 
anniversary of the discovery of radio by Marconi, Sep- 
tember 28 to October 5, 1947, in Rome, Italy. The event 
is being organized by the Italian National Council of 
Research. It aims at providing a complete picture of the 
progress of radio in the world from the scientific, technical, 
and industrial point of view, and giving world radio experts 
the opportunity of meeting and studying plans for a future 
international collaboration. 
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The Bulletin of Mathematical Biophysics 


Following is the table of contents for The Bulletin of 
Mathematical Biophysics for June 1947: 
Mathematical Theory of Motivation Interactions of Two Individuals: 
II. ANATOL RAPOPORT 


A Theory of Membrane Permeability: III. The Effect of Hydrostatic 
Pressure. INGRAM BLOocH 


The Mechanism of the Middle Ear: Il. The Drum. Martinus H. M. 
ESSER 


A Mathematical Description of Metabolizing Systems: I]. HerMAN 
BRANSON 


Outline of a Matrix Calculus for Neural Nets: II. Il. D. LANDAnI 
The University of 
Volume 9, Number 2. 


Chicago Press, Chicago, Illinois, 


Meeting of the Division of Solid State Physics 


\ Division of Solid State Physics in the American 
Physical Society, recently authorized by the Council, is in 
the process of formation and a corresponding announce- 
ment is being sent out to the members of the society. The 
organizing committee of the division is planning for a 
meeting of the division in conjunction with the Montreal 
meeting of the American Physical Society June 19 to 21, 
1947. In view of the fact that the American Society for 
X-Ray and Electron Diffraction is having a meeting near 
Montreal starting Monday, June 23, the meeting of the 
division will take place on Saturday, June 21, so as to 
make it convenient for those interested in solid state to 
attend both meetings. 

The program will be composed of four invited papers: 
“Diffraction of Neutrons in Crystals,"’ F. Seirz, 

“Theory of Crystal Rectifiers with Application to Copper Oxide,” 

J. BARDEEN, 

“The Relaxation Spectrum of Metals,’ C. ZENER, 
“Nuclear Magnetic Resonance Absorption in Solids and Liquids,” 

E. M. PURCELL, 


and contributed papers in the field of solids. It is suggested 
that those who have papers in the field of the division 
would present them at the June meeting in Montreal. 
The Organizing Committee of the Division of Solid State: 
T. A. Read S. Siegel 
F. Seitz R. Smoluchowski 
W. Shockley C. Zener 





Calendar of Meetings 


Institute of Aeronautical Sciences, Detroit, Michigan 
American Physical Society, New England Section, 
Massachusetts 


Amherst, 


Society of Automotive Engineers, French Lick, Indiana 

American Society of Refrigerating Engineers, Los Angeles, 
California 

American Institute of Electrical Engineers, Montreal, Quebec, 
Canada 

Symposium on Molecular Structure and Microscopy, Ohio 
State University, Columbus, Ohio 


12-14 Canadian Association of Professional Physicists, University of 
Western Ontario, London, Canada 

16-19 American Society of Mechanical Engineers, Chicago, Illinois 

16 American Society for Testing Materials, Atlantic City, New 
Jersey 

18-2 American Association of Physics Teachers, Regional Meeting, 
University of Minnesota, Minneapolis, Minnesota 

i8 Society for the Promotion of Engineering Education, Minne- 
apolis, Minnesota = 

23-25 American Society for X-Ray and Electron Diffraction, Mon- 
treal, Quebec, Canada 


July 
11-12 
16-19 


American Physical Society, Stanford University, California 
American Society of Civil Engineers, Duluth, Minnesota 


August 


21-23 Society of Automotive Engineers, Los Angeles, California 
(Transportation and Maintenance Meeting) 
American Institute of Electrical Engineers, San Diego, Cali- 


fornia 


26-29 





New Booklets 








Dow Corning Corporation, Midland, Michigan, has 
announced a new technical pamphlet, DC Antifoam A, 
which discusses a silicone compound recently developed 


for killing foam in aqueous systems. 4 pages. 


The December 1946 issue of The Frontier, published by 
the Armour Research Foundation of Illinois Institute of 
Technology, Chicago 16, Illinois, includes the following 
articles: “A Unique Torquemeter,” 
Analytical Methods,” 


Refractory Electrical Insulating Compositions.” 


“Classification of 
“Processing and Use of Zircon for 
20 pages. 


Nickel Steel Topics for February 1947, published by 
The International Nickel Company, Inc., 67 Wall Street, 
New York 5, New York, article on the 
“refrigerator car of the future,”’ built largely of stainless 
and low nickel alloy steels. 12 pages. 


features an 


Fish-Schurman Corporation, 230 East 45th Street, New 
York 17, New York, makers of raw optical, technical, and 
scientific glass and optical parts and specialties, recently 
MI-318 on FS Multi-Layer Interference 
Films; Dichroic and Achromatic Beam Splitters. 2 pages. 


issued bulletin 


The Aerovox Research Worker, monthly house organ of 
the Aerovox Corporation, New Bedford, Massachusetts, 
in its 


December 1946 edition discusses ‘‘The Rectifier- 


Amplifier V. T. Voltmeter.” 4 pages. 


Pittsburgh Plate Glass Company, Columbia Chemical 
Division, Fifth Avenue at Bellefield, Pittsburgh 13, Penn- 
sylvania, has published Form A-100, a new Columbia 
Caustic Soda Book, 72 pages, printed in three colors. 
Many photographs. Write for a copy on your company 
letterhead. 
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